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1. DESCRIPTION OF RESEARCH AND RESULTS 



The research goal of this project may be summarized into the 
following three parts 

a) Waves forced by condensation heating in the troposphere 

We have studies the problem of the vertical structure of trop- 
ical waves generated by cumulus heating for two cases: l)heating 

is controlled by large-scale low-level convergence (CISK model) 
and 2) heating receives no feedback from large-scale motion 
(forced model) . In the first case the internal wave-CISK system 
is treated and upper and lower bounds for the vertical wavelength 
of the unstable waves under normal heating conditions are esta- 
blished through analysis of the frequency equation. The lower 
bound excludes the possibility of excitation or maintenance of 
short vertical wavelengths relative to the vertical scale of heat- 
ing. The vertical structure of the most unstable waves is also 
computed. These results were reported in Chang (1976a) . In the 
second case the problem of scale-selection of Kelvin waves in the 
stratosphere is analyzed using an inviscid linear model. The re- 
sults modify Holton's (1973) theory in that 1) the forcing is 
most efficient for the longest zonal wavelength even if the heat 
sources are distributed randomly, and 2) the most favored verti- 
cal wavelength of the excited waves is about twice the vertical 
scale of heating. These results were reported in Chang (1976b). 

We have also treated the linear model with friction added which 
is called for by several budget studies of cumulus -related tropi- 
cal weather systems, including the monsoon (Holton and Colton, 
1972), the cloud clusters (Williams and Gray, 1973), and the 
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easterly waves Recker, 19 71; Wallace, 19 71, Reed and 

Johnson, 1974; Chu, 1976]. This friction is believed to be due 
to vertical cumulus transports. Our analysis shows that this 
damping has a strong influence at low frequencies on the forced 
equatorial waves and results in two types of dispersive relation- 
ships. The first type is characteristic of the regular internal 
gravity waves which have fast phase speeds and weak vertical atten- 
uation. The second type is dominated by the viscous damping time 
scale and has slow phase speed and strong vertical trapping. This 
theory predicts that the stratospheric oscillations may be identi- 
fied with the first type and the tropospheric oscillations with 
the second type. In the case of Kelvin waves the results can be 
used to explain consistently both the observed stratospheric Kelvin 
waves and the planetary-scale Kelvin-like oscillations in the 
troposphere including the 40-50 day oscillation and the winter mon- 
soon and Walker circulations. These results are described in 
Chang Cl977a). 

b) The local barotropic instability produced by the spatially 
varying basic flow . 

Here we have studied the behavior of waves superimposed upon 
a barotropically unstable mean wind field, which varies in x and 
y. This mean wind field roughly simulates the 200 mb easterly jet 
south of the Tibetan high during the northern hemisphere summer. 

The barotropic vorticity equation is linearized and a linear damp- 
ing term is added. The following mean u component is used 

uC.x,y] = -UCx) sech^ Cy/d(x) ) - , 
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where dCx) is a measure of the width of the jet. The quantity 
U(x) and u are chosen in such a way that the mean flow is non- 
divergent. A rectangular domain is employed with a time periodic 
boundary contition on the east and a radiation condition on the 
west, which simulates the propagation of small amplitude waves 
through the easterly jet region. The vorticity equation is solved 
with the use of finite differences, and when the boundary condi- 
tions are properly adjusted the waves move smoothly across the re- 
gion and out the western boundary. After a certain time interval 
of numerical integration, the solution becomes periodic everywhere, 
with the frequency that is specified on the eastern boundary. As 
the waves move through the domain they grow or decay spatially in 
reaction to the local basic flow, even though at each point the 
variation is purely harmonic. We have developed a simple mechan- 
istic analytical model which describes the main features of the 
spatial variation. The simple model uses growth rates and phase 
velocities which are computed from a parallel flow model with the 
wavelengths taken from the full numerical model. These wave- 
lengths increase to a miximum at the point where the jet is most 
unstable, and then decrease; the phase velocity has a similar be- 
havior. The maximum amplitude in the numerical solution occurs 
downstream from the most unstable part of the jet, at the point 
where the growth rate in the parallel flow model becomes zero. 

Tire above results are grossly similar to that expected from the 
parallel flow theory of barotropic instability. However, in the 
unstable region the resultant structure of the waves causes a spa- 
tial growtli rate greater than predicted by the local growth rates 



7 



computed with a parallel flow model. In the stable region, the 
structure leads to a strong dynamic damping which is apparently 
due to the presence of the continuous spectrum modes CCase, 1960; 
Pedlosky, 1964). When a uniform advective velocity is added to a 
variable mean flow, the difference between the behavior of the 
computed waves and that implied by the parallel flow theory is 
somewhat reduced. However, in this case a stronger zonal asymmetry 
in the spatial growth rate with respect to the jet maximum occurs 
as a result of slower adjustment of the wave structure to the local 
stability conditions. Chang and Williams (1976) presented the 
analytic model and discussed some other analytical aspects of the 
problem. The numerical model results and comparisons with the 
analytic model will be reported by Tupaz , Williams and Chang (1978). 

c) The numerical modeling of the wave-basic flow interacting 
system . 

In this part we adopted a simple numerical model formulated 
by Monaco and Williams (1975) to simulate the planetary-scale mon- 
soon circulation and the interaction with synoptic waves. In 
order to concentrate on the vorticity dynamics, the heating is 
specified with a horizontal distribution given by the time-mean 
200-mb divergence of the 1967 summer analyzed by Krishnamurti and 
Rogers (1970). The domain of the model extends from 18°S to 
46°N and there are three o-levels in the vertical. Topography 
is not included although its thermal effect is implied in the 
heating function. Two preliminary experiments have been carried 
out to study the vorticity budget of the planetary scales and the 
non-linear barotropic energy conversions. In one experiment the 
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heating is fixed in time while in the other it fluctuates with a 
period of 10 days. The time-mean fields with moderate damping in 
both experiments show a significant westward phase shift of the 
Tibetan high, which suggests that the correlation of the fluctua- 
tions of the divergence and vorticity as suggested by Holton and 
Colton (1972) is probably not sufficient for the required planetary- 
scale vorticity sink. However, the integration period is only 
three times the fluctuation period in the second experiment and 
the amplitude of fluctuation is about half that observed. So the 
full effect of the fluctuation may not have been properly repre- 
sented in the experiment. A longer integration with larger fluc- 
tuation is probably necessary to assess fully the effect of the 
fluctuating heating. On the other hand, the positions of the 
simulated tropical upper tropospheric troughs in our model agree 
quite well with those observed. This may suggest that damping due 
to cumulus transport is more likely to account for the vorticity 
sink in the Tibetan high area. V/e also found singificant transient 
wave activities at the upper levels in both experiments and they 
appear to he drawing energy from the planetary scale barotropi- 
cally, in accordance with the observation by Krishnamurti Cl^^l) . 

It is possible that they play an important role in reducing the 
amplitudes of the long waves which agree quite well with those 
observed despite the phase problem. This point is borne out by 
comparison with the excessively large amplitude found in Holton 
and Colton's linear model when strong damping is excluded. In 
addition, these short waves are most active in the fluctuating 
heating case, suggesting that local barotropic instability at 



9 



upper levels may be more important than that implied by the time- 
mean wind. These results have been reported by Chang and V/illiams 
(1976b) and Chang and Pentimonti (1977). 

Some of the results concerning the planetary-scale monsoon, 
including the interpretation of the stationary equatorial winter 
monsoon circulation as a viscous, thermally-forced Kelvin wave 
(part a) , the preliminary result of the study of barotropic in- 
stability of the upper- tropospheric summer monsoon basic flow 
(part b) , and the preliminary result of the planetary scale vor- 
ticity budget using forced numerical model (part c) are included 
in a review paper by Chang (1977b). The research on these prob- 
lems is being continued under NSF Grant ATM77- 14821 which began 
on 15 September 1977. 



10 



REFERENCES 



Case, K.M. , 1960: Stability of inviscid plane Couette flow. Phys 

of Fluids , 3, 143-157. 

Chang, C.-P., 1976a: Vertical structure of tropical waves main- 

tained by internally-induced cumulus heating. J. Atmos. Sci., 
32, 729-739. 

Chang, C.-P., 1976b: Forcing of s tratopsheric Kelvin waves by 

tropospheric heat sources. J . Atmos . Sci . , 33 , 740-744. 

Chang, C.-P., 1977a: Viscous internal gravity waves and low- 

frequency oscillations in the tropics. Submitted to J. Atmos. 
Sci. 



Chang, C.-P., 1977b: Some theoretical problems of the planetary- 

scale monsoons. Pure and Appl. Geophys . , 115 , 1089-1109. 

Chang, C.-P. and R.T. Williams, 1976a: Barotropic instability of 

a spatially- varying mean flow. Conference on atmospheric 
waves and stability, Seattle, March-April 1976. Bull. Amer . 
Meteor . Soc . , 5 7 , 101. 

Chang, C.-P. and R.T. Williams, 1976b: Study of planetary and 

synoptic-scale waves of the northern summer monsoon. CARP 
Study Conference on tropical numerical modeling, Exeter, 
England, April 1976. 

Chu, J.-H., 1976: Vorticity in maritime cumulus clouds and its 

effects on the large-scale budget of vorticity in the tropics. 
Ph.D. thesis, UCLA, 123 pp . 

Holton, J.R., 1973: On the frequency distribution of atmospheric 

Kelvin waves. J . Atmos . Sci . , 30 , 499-501. 

Holton, J.R. and D.E. Colton, 1972; A diagnostic study of the 

vorticity balance of 200 mb in the tropics during the northern 
summer. J. Atmos . Sci . , 29 , 1124-1128. 

Krishnamurti , T.N., 1971: Tropical east-west circulations during 

the northern summer. J . Atmos . Sci . , 2 8 , 1342-1347. 

Krishnamurti, T.N. and 3. Rogers, 1970: 200 millibar wind field 

June, July and August, 1967. Pv.ept 70-2, Dept, of Meteorology, 
Florida State University, Tallahassee. 

Monaco, A.V. and R.T. Williams, 1975: An atmospheric global pre- 

diction model using a modified Arakawa differencing scheme. 
Naval Postgraduate School Technical Rept. NPS- 51Wu75041 , 86 pp 

Pedlosky, J. , 1964: An initial value problem in the theory of 

barotropic instability. Tellus , 16 , 12-17. 



11 



Reed, R. J. and R. H. Johnson, 1974: The vorticity budget of sy- 

notpic-scale wave disturbances in the tropical western Pacific, 
J. Atmos. Sci., 31, 1784-1790. 

Reed, R.J. and E.E. Recker, 1971: Structure and properties of sy- 

noptic-scale wave disturbances in the equatorial western Pacific 
J. Atmos. Sci ., 2 ^, 1117-1133. 

Tupaz , J.B. , R.T. Williams and C.-P. Chang, 1978: A numerical 

study of barotropic instability in a zonally varying easterly 
jet. J . Atmos . Sci . , 35 , C.in press). 

Wallace, J.M. , 1971: Spectral studies of tropospheric wave dis- 

turbances in the tropical western Pacific. Rev. Geophys . 

Space Phys . , 9_, 557-612 . 

Williams, K.T. and W.M. Gray, 1973: Statistical analysis of 

satellite-observed trade wind cloud clusters in the western 
North Pacific. Tellus , 25, 313-336. 



LIST OF PUBLICATIONS 



Chang, C.-P., 1976a: Vertical structure of tropical waves main- 

tained by internally- induced cumulus heating. J. Atmos. Sci., 
33, 729-739. 

Chang, C.-P., 1976b: Forcing of stratospheric Kelvin waves by 

tropospheric heat sources. J . Atmos . Sci . , 33 , 740-744. 

Chang, C.-P., 1976c: Comments on "Instability theory of large- 

scale disturbances in the tropics." J. Atmos. Sci., 33, 
1667-1668. 

Chang, C.-P., 1977a: Viscous internal gravity waves and low- 

frequency oscillations in the tropics. Submitted to J. Atmos . 
Sci . 



Chang, C.-P., 1977b: Some theoretical problems of the planetary- 

scale monsoons. Pure and Appl. Geophys ., 115 , 1089-1109. 

Chang, C.-P. and R.J. Pentimonti, 1977: A numerical study of 

time-mean northern summer monsoon with steady and fluctuating 
heating. Proc. Internatl. Symp . on Monsoons, New Delhi, India. 

Tupaz, J.B., R.T. Williams and C.-P. Chang, 1977: A numerical 

study of the locally unstable barotropic easterly jet. Proc. 
Internatl. Symp. on monsoons. New Delhi, India. 



12 



THESES 



Pentimonti , R.J.: Study of monsoon circulation with steady and 
fluctuating heating. M.S. thesis. December 1976. 

Tupaz , J.B.: A numerical study of barotropic instability of a 

zonally varying easterly jet. Ph.D. thesis, June 1977. 



SCIENTIFIC COLLABORATORS 

R.J. Pentimonti, graduate student 

J.B. Tupaz, graduate student 

L.C. Chou, graduate research associate 



COMMENTS 



This research is being continued under NSF Grant ATM77-14821. 



13 



Reprinted from Journal of the Atmospheric Sciences, Vol. 33, No. S, May 1976 
American Meteorological Society 
Printed in U. S. A. 



Vertical Structure of Tropical Waves Maintained by 
Internally-Induced Cumulus Heating 

Chih-Pei Chang 

Department of Meteor ology^ Naval Postgraduate School ^ Monterey^ Calif. 93940 
(Manuscript received 25 November 1975, in revised form 22 January 1976) 

ABSTRACT 

Solutions to the wave-CISK (conditional instability of the second kind with cumulus heating being in- 
duced by low-level internal wave convergence) system are obtained to study the vertical structure of margin- 
ally unstable waves. A diabatic heating profile is specified that resembles those observed and those theoret- 
ically derived from simple parameterization schemes. Upper and lower bounds for the vertical wavelength 
of the unstable waves under normal heating conditions are established through analysis of the frequency 
(stability) equation. The lower bound excludes the possibility of excitation or maintenance of short vertical 
wavelengths (relative to the vertical scale of heating) by wave-CISK. The calculated growth rates indicate 
that this result is basically insensitive to the vertical heating profile. The vertical structure of the most 
unstable waves is also computed and the possible roles played by CISK in large-scale tropical waves are 
discussed in light of these results. 



1. Introduction 

The possibility that large-scale waves with short 
vertical wavelengths may be excited in the tropical 
atmosphere has been raised by Lindzen (1967), who 
showed that tfie equatorially trapped, internal gravity 
and internal Rossby waves usually have a small 
positive equivalent depth and therefore a short ver- 
tical wavelength. Holton (1969, 1972a) later pointed 
out that such short vertical wavelengths, if they exist 
in the troposphere, would cause great difficulties for 
numerical weather prediction in the equatorial region 
because of the stringent requirement for high vertical 
resolution. Whether short vertical-scale waves actu- 
ally have a significant presence in the tropical tropo- 
sphere or not is unclear due to inadequate observations. 
Studies based on conventional radiosonde data (e.g., 
Wallace, 1971; and others) tend to suggest that most 
of the variance of the waves in the tropical tropo- 
sphere is associated with vertical scales comparable 
to the scale height. However, Madden and Zipser’s 
(1970) analysis of high- resolution rawinsonde data 
obtained during the Line Island Experiment showed 
that a quite short vertical wavelength (^^3 km) 
existed over the central equatorial Pacific. Further 
observational studies thus seem necessary to resolve 
this problem. 

Condensation heating due to deep cumulus con- 
vection has been recognized as a primary energy 
source for large-scale tropical motions. Two possible 
mechanisms which may be responsible for the excita- 
tion of large-scale tropical waves by heating have 



been investigated. The first is a one-way forcing by 
the cumulus scale which receives no feedback from 
the large scale. Numerical study by Holton (1971, 
1972a, 1973) and analytical work by Chang (1976) 
have shown that this mechanism excites tropical 
waves that have vertical scales comparable to those 
of the forcing, which are usually 10-14 km. The 
problem of short vertical wavelengths thus appears 
unimportant for forced waves. The second possible 
mechanism is the conditional instability of the second 
kind (CISK) which involves a feedback to cumulus 
convection by low-level convergence of large-scale 
motion. This low-level convergence may be due to 
frictionally induced cross-isobaric flow, as in the case 
of t>T>hoons (which shall be called ‘‘Ekman-CISK”), 
or due to the structure of the internal waves (which 
shall be called ^Vave-CISK’O* The Ekman-CISK has 
been shown by Chang (1971) and Chang and Piwowar 
(1974) to possess no preferred horizontal scale for 
synoptic wave disturbances. Thus Ekman-CISK may 
not explain the generation of synoptic-scale waves, 
but it may be quite important in maintaining the 
waves once they are initiated. The wave-CISK mecha- 
nism has been investigated by Hayashi (1970), Lindzen 
(1974) and Lindzen et al. (1975). Although Lindzen 
(1974) and Lindzen et al. (1975) found no direct 
horizontal scale selection, the most unstable wave- 
CISK mode has a vertical scale seemingly comparable 
to the depth of the sub-cloud layer or the low-level 
convergent moist layer, depending on the particular 
model aspects. A short vertical scale may therefore 
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be important for the excited waves. Using the cal- 
culated vertical scale they also speculated that, due 
to the Doppler-shifting effect of the variable basic 
flow and the small growth rates, a period of ^^5 days 
associated with a wavenumber zero, mixed Rossby- 
gravity mode is likely to dominate the tropical spectra. 
This periodicity thus provides a basic frequency for 
tropical systems. Depending on the local basic flow, 
it may force resonantly higher wavenumbers that may 
be relevant to the observed waves. 

Since the growth rates associated with wave-CISK 
are usually very small, it may be expected that the 
CISK waves have properties similar to neutrally 
forced waves. The difference in the selection of vertical 
scales between the forced models and the wave-CISK 
models thus appears to be worth further investigation. 
In addition to its importance for tropical numerical 
weather prediction, the vertical scale also determines 
other properties of the excited internal waves. The 
purpose of this paper is to examine the growth rates 
associated with W'ave-CISK, particularly with regard 
to the vertical scale selection mechanism. 



vertical structure equations by assuming that 
w = '£, Vniy)'w„{z) exp[s/(2//)], 

n 

where the density factor exp [z/ (2/7)] is separated 
for convenience. The meridional structure equation is 

dWn /k0 5y\ 

+( k^+ )Fn=0, (6) 

dy^ \ oi ghn ghj 

where g is the gravitational constant and hn, the 
equivalent depth, is the separation constant. Matsuno 
(1966) and Lindzen (1967) have shown that the solu- 
tions to (6) which satisfy the boundary condition 

Fn“»0 as ly|— (7) 

and lead to the frequency equation 

fk0 02^\{ghn)'^ 

{ F+— = 2«+l, n=~l,0,l,2, (8) 

\ w ghJ 0 



2. Basic equations 

Neglecting the shear of the mean zonal wind, the 
linearized zonal and meridional momentum, hydro- 
static, thermodynamic energy and continuity equa- 
tions for a single zonal wavenumber k on an equatorial 
beta-plane may be written 



are 



Yn{y) = 






1 ighn)'^ 

(Ji 



2[l+^(g*n)5] 






Xexp(-?V2), 



ioiU —jdyv = ~ ik<py 


(1) 


d<p 

i(j3V-\-^yu— , 

dy 


(2) 


d4> Rt 
dz H ' 


(3) 


iut+wT=—, 

Cp 


(4) 


dv d 

iM-\ — (e"^^^iy) = 0, 


(5) 



dy dz 

where zl, v, w, T, $ and Q are the perturbation zonal 
velocity, meridional velocity, vertical velocity, tem- 
perature, geopotential and diabatic heating rate, 
respectively; w is the Doppler-shifted frequency, 
H a constant scale height, T the static stability, 
Cp the specific heat at constant pressure, R the gas 
constant, 0 the meridional gradient of the vertical 
component of earth’s vorticity, y the meridional co- 
ordinate, and s= —// In {p/po) is the vertical coordi- 
nate with p the pressure and /?o a reference pressure. 

Eqs. (l)-(5) may be combined into a single equa- 
tion in w that may be separated into meridional and 



where and the Hermite polynomials 

Hn{i) have values only for The meridional 

velocity solution is ”t^n(y) °^//n(^) exp (~^/2). Lamb 
(1973) has solved (6) using the more general form 
of the confluent hypergeometric functions without the 
restriction of (7), In such a case, the lateral boundary 
conditions are given by the choice of vanishing me- 
ridional velocity at some finite (y| where the con- 
fluent hypergeometric functions have a node. The 
Hermite polynomials are a sub-set of the confluent 
hypergeometric functions. 

The vertical structure equation which is most rele- 
vant in the present calculation is 

dhjOn Qn 

h ) 

dz^ ghn 

where 

S 1 

Xn^ = , 

ghn 4/72 



R 




R ^ 

Qn= QnYn~'- eXp[-s/ {2H)~\. 



(9) 

( 10 ) 
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The parameter Xn is a measure of the vertical wave- 
number and is generally complex for unstable waves. 
In the case of forced modes a; and k are known pa- 
rameters of the forcing function so that hn and Xn are 
obtained from (8) and (10) respectively. For CISK 
modes, on the other hand, Xn is determined as an 
eigenvalue of (9) and oj is then found from (8) and (10). 
In the present calculation, the heating function Q' is 
assumed to possess a white-noise distribution for all 
Hermite modes (all values of n) so that the calcula- 
tion of h will be independent of n. The subscript n 
will therefore be dropped in the subsequent discus- 
sions. It is well known that the vertical structure 
equation (9) can also be derived by considering two- 
dimensional wave motion only. 

As in the forced model by Chang (1976) the fol- 
lowing boundary conditions are used to solve (9) : 



w 



0, at z=0 

at z — zt 



(lla) 

(llb) 



where 



Cl — 



Here Zt is the height of tropopause, Ci, Co are con- 
stants, and the condition (lib) results from the re- 
quirement that latent heating vanishes at The 
parameter r is a reflection coeflicient which, in the 
absence of vertical wind shear, is given by the speci- 
fication of the static stability distribution. If the 
heating is assumed to vanish below the cloud base (zc), 
the solution to (9) may be written in the form 



w{z) — / sinXz — dz, z'^Zt, (12a) 

X(l-fr) Ac gh 



sinXz Q' 

zi'(z) = / — dz 

X(l+r)j, gh X(l+r) 



r Q' 

X / sinXz — dzy 
J gh 



Zt>z>Zcy (12b) 



sinXz Q' 

w{z)=^ / — dz. 

X (!+'') Ac gh 



Zc^z^O. (12c) 



The solution procedure using the Green’s function 
technique is the same as in Lindzen (1974). 



3. The heating function 

The parameterization of the effect of condensation 
heating due to cumulus convection is a major problem 
in tropical modeling. Important contributions on this 
problem have been made by the sophisticated schemes 
developed by Arakawa and Schubert (1974), Ooyama 
(1971), and others. On the other hand, in simplified 



analytical models, simpler schemes must be used to 
keep the mathematical analysis tractable. Usually in 
these models the vertical distribution of the heating 
is specified as a given function that is invariant in 
time. The major criticism of this type of parameteriza- 
tion is that the results are sensitive to the specified 
heating profile such that one can obtain the desired 
result by simply tuning the heating profile. Another 
criticism is that the inherent time-variation of the 
heating function due to the complicated scale inter- 
action cannot be properly included. However, recent 
observational studies, using both the spectral analysis 
technique and the composite technique, have produced 
a consistent picture of the vertical heating profile due 
to latent heat release over the tropical oceanic areas. 
The profiles observed by Nitta (1970), Wallace (1971), 
Reed and Recker (1971) and Williams and Gray (1973) 
are plotted in Fig. 1. It can be seen that in each case 
a maximum occurs in the middle troposphere between 
the 6 and 9 km levels. On the other hand, observa- 
tional analysis of the divergence and vertical velocity 
by Schubert and Reed (1975) using recent GATE 
data has yielded a result that implies a vertical heating 
profile different from those shown in Fig. 1. They 
found that the maximum vertical velocity, which 
usually coincides wfith the maximum heating, based 
on scaling arguments (Holton, 1972b), occurs at a 
much lower level. How^ever, analyses of the same data 
set by Petrossiants ei aL (1975) and Antsipovich et aL 
(1975) of the USSR Hydrometeorological Center have 
resulted in profiles which are in close agreement wfith 
Fig. 1. In any case, the profiles included in Fig. 1 are 




Fig. 1. Vertical heating functions deduced from observations and 
from simple theoretical parameterization schemes. 
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based on a much larger sample of data and may be 
viewed as an indication that the tropical atmosphere, 
on a statistical basis, tends to maintain a similar 
vertical profile for the large-scale heating function. 

Results of a theoretical calculation by Kuo (1965) 
using his parameterization scheme are also included 
in Fig. 1. Kuo’s parameterization was originally based 
on the assumption of horizontal mixing between cloud 
air and environmental air. Recently Kuo (1974) has 
shown that his formulation is equivalent to the con- 
sideration that heating is due to the adiabatic com- 
pression of the descending environment. Kuo’s scheme 
has been used in many tropical numerical models 
with a degree of success not inferior to any other 
parameterization scheme. It is therefore interesting to 
compare his heating profiles with those observed. In 
Fig. 1, the Kuo-1 curve is for the case of no entrain- 
ment and the Kuo-2 curve is for the case of moderate 
entrainment. It is clear that the latter curve, having 
a maximum near 7 km, agrees quite well with those 
deduced from observations. The average heating profile 
obtained from the GFDL general circulation model 
as calculated by Hayashi (1973) is also included, 
because the GFDL model uses a different parame- 
terization scheme: the moist convective adjustment. 
Hayashi decomposed the heating function into com- 
ponents of various equatorial beta-plane modes. Among 
them, the profile for the Kelvin waves (>i= — 1) at 
the equator is one of the curves least resembling 
those observed. This profile is also plotted in Fig. 1. 
It is seen that even for this highly simplified parame- 
terization scheme the maximum heating still occurs 
near 6 km, within the range of the other curves. 




Fig. 2. The specified heating function Q'{z). 



Based on the above discussion the heating functior le 
for the wave-CISK parameterization is specified as 



( 2 ) = \mNStWi,e°'^' sin tts ', Zt'^z^Zc 






z>zt or z<Zc 



(13: 



where 



S — 2c 



S =- 



Az 



and ^z — Zt^Zc is a measure of the vertical scale oi' 
the heating. The coefficient m specifies the strength 
of the heating and the factor \ is needed becausF^ 
the amplitude of the Fourier component is one-hal|e 
of the heating maximum, as only positive conden 
sation heating is permitted. The heating is proporliti 
tional to the moisture convergence in the moist layel^fc 
which is represented by the vertical velocity Wh at th^ie 
top of the moist layer. The tropospheric value of th^iaj 
static stability St is included in the proportionalit} 
constant. The parameter a is used to vary the maxi-’J^ 
mum heating level to test the sensitivity of the model’um 
The coefficient N is a normalization factor so thatite 
when other conditions are equal the total amount 
heat release weighted by density in a column wouldfu 
remain the same with different values of a. Thus, 



.V=- 



f exp[-s, 
Jo 



,li2H)2dz' 



/ exp^az' — 

J 0 



z/(2H)'] sinirz'dz' 



2/7 { 1 - exp[ - Az/ (2//) ] } { [a - Ac/ (2H)J+it^} 






T, 



:pr 



k 



7T { exp[a — Az/(2//)]-fl} 



Profiles of (13) multiplied by the density factor 
exp ( 2 / 2 //) are plotted in Fig. 2 for various values 
of a. The range of — with maximum heating 

between 7 and 9 km, appears to be representative of 
the observed and theoretical profiles shown in Fig. 1, 



.na 



4. Analysis of the stability equation 



If the mixed layer with top at Zb (<Zc) is assumed 
to be the moisture layer which provides the con-; 
vergence for CISK, the most relevant solution in 
studying the instability is (12c) evaluated at Zb. The 
solution is 






111 ;; 

'>eni 



sinXzi mwi St irAz 



’4: 

M 



WI,= 



\+r gk 2 



q- 



X 












where m = mN and q=a — i\^z. Equating the coeffi- 
cients of Wb on both sides and using (10) leads to 



\1 
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e stability equation 



[\)^-mZbSz 

2 



( 



AIPJ 



X- 



(1+0 (?■++) 



-1=0, 

(14) 



gives two equations : 

Re(^')=0, 
lm(^') = 0. 

For <2 = 0, (18) may be expanded to give 

Re('0 = fx.2~Xr+— Vexp(X.3,)cosX.s, 

\ 

+exp(X,Sf)cosXr2 c]+2Xr/\,Cexp (Xi2^)sinXr 



(18) 



lere we have assumed that (sin X26)/X~26. It can 
shown that g-+ 7 r" = 0 is not a solution to (14). 
p (14) must be solved for the complex eigenvalue X 
ing either numerical or graphical methods. However, 
fore doing so some valuable information regarding 
e range of vertical scale of the unstable solutions 
ly be obtained by analyzing (14). 

The complex eigenvalue X may be written as X=Xr 
fXj, where the real part Xr is the vertical wave- 
imber while the imaginary part Xi gives the growth 
te of the waves. If the positive Xr is chosen, the C 2 
mponent in (11b) represents the downward phase 
upward energy propagation and the Ci component 
presents the upward phase or downward energy 
opagation. It then follows from (10) and (8) that 
X, is negative, the imaginary part of the frequency cj 
11 be negative and the waves will grow exponentially 
time. Based on observations, and justified a pos- 
ion\ we shall consider only the case of marginalh' 
(Stable waves, or very small growth rates. Thus the 
nditions applied to (14) may be written 

X.= ~iX*|, (15) 

X,2«X,2. (16) 



. addition, we shall assume 



1 

Xr'» , 



(17) 



^ich is valid for all practical purposes since the 
ale height H of the atmosphere is ^7 km and (17) 
ily places a restriction of ^^88 km for the maximum 
Ttical wavelength. 

Eq. (14) will now be analyzed b}’ assuming that 
=0 and that the stratosphere {z>Zt) and troposphere 
^ 2 <) have the same constant static stability St. 
[ this case r = 0 and since g^+7r-?^0 the stability 
[uation may be re-written as 



'(X) =4^(X) (^=+7T=) =-m26A/x2+— Ve^+1) 

2 \ 4//V 

X exp ( — iXtZc) = 0. 

iparating 4^'(X) into its real and imaginary parts 



2[7t2--A22(X.2~X,2)] 

+exp(XxsJsinX,Sc]d ^ = 0. 

TT'mZb^Z 

Simplifying this equation with the conditions (16) 
and (17) and using (15) lead to 

Re('F') 

- (x;^ cosXrCr — 2Xr!X,!sinX,Si)exp(— IXejsc) 

I 

+ (Xr^ COsXr^c — 2Xr! Xe | sinXr2c)eXp(- ' X, | 2^) 



II 

2(7T— As2X,2) 

+ = 0. (19) 

TT7nZbSz 



III 



Since | cos Xr 2 < ; , sin Xr 2 ^ ' , cos XrZc and sin XrZc 
are all ^1 and exp(~ * Xj ' 2^) < 1, from (16) we have, 
allowing for small error under extreme cases,^ 



IT^X.^ 

|II|^X.2 



( 20 ) 



Thus for (19) to be satisfied we may write 

Max[Re(4'')] = Max(I)+Max(II)+III>0, 

where Max denotes the maximum value using the 
range (20) : 

Max(I) = Max(II)=X.l 



Thus we have 



which yields 



if 



27t 2AsX.2 

2X.2+~ —>0, 

mzbSz TrrnZb 



Xr-< 

Az{Sz — 7rmZb) 

Az>wmzb. 



( 21 ) 

( 22 ) 



^ The precise maximum for terms I and II is ( | (l-|-l€-)cos 
Qtan“*(~2€) I ] exp(— |X, |cc) iXr“, where e=|X,/Xrl. Inequality 
(20) is true when € — 0. For e = 0.1 which is our assumption the 
maximum is [1.02 exp(— IXi |cc)]Xr^. Obviously if ^\{\ becomes 
large the exponential factor would become much less than 1. 
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k=1 ks2 k=3 k>4 k-5 (<^6 




Fig. 3. Values of Re('I^) and Im('I') as a function of Xr[km and Xi[km The crosses in the ReW 
diagram are solutions for which ReW = lm(^) = 0. 



Choosing the minimum values (Min) in (20), we have 

Min[Re(^')] = Min(I)+Min(II)+IIKO, 
or 

27t 2AsX,2 

-2.V+- 1— <0, 

jnzhi^z Ttmzb 

giving 

^2 

Xr'> . (23) 

Az{Az+miZb) 



so the vertical wavelength, L — lw/Xn is limited by 
the following bounds : 

13 km<L<40 km. 

Thus waves with a vertical wavelength shorter than 
13 km will not be excited by wave-CISK. If such 
waves are first excited by another mechanism, the 
above result indicates that they cannot be mahitained 
by wave-CISK. 



So (21) and (23) give the upper and lower bounds 
for Xr. It is interesting to note that Xt-^tt/As in the 
limit of /7?— >0, corresponding to a vertical wavelength 
of twice the heating scale, which is the most efficiently 
excited wavelength for forced waves (Chang, 1976). 
However, Xr = 7 r/A 2 cannot be a solution for any 
finite fh because it does not satisfy (14). 

Now we shall use the values 24=15 km, 2^=1 km, 
A2=14 km, 26 = 0.4 km, m = 5.8 which leads to rn = 9. 
The value of m is calculated by fitting (13) to Reed 
and Recker^s (1971) observed heating and vertical 
velocity profiles for the wave trough, for which an 
estimate of 

m 26 — 2.32 km, (24) 

can be made assuming that 5*4=1.22X10“^ s“^ based 
on a tropospheric value of F = 3°C km“L Eq. (24) is 
a good approximation as long as 2&<4 km, so the 
particular choice of 26 is not crucial. From Reed and 
Recker^s data z£;6=^0.29 cm s“^ so that the maximum 
heating is ^6.8°C day”h Substituting these values 
into (21)-(23) we obtain 

7T 7T 

= 0.512> |Xri > = 0.167 km“\ 

6.14 18.8 



5. Eigenvalues and eigenfunctions 

a. Growth rates 

The stability equation (14) will now be solved using 
the graphical method with the accuracy of the solu- 
tions Improved by a numerical iteration technique. 
The graphical method is illustrated in Fig. 3, where 
the real and imaginary parts of ^(X) given by (14) 
with r = 0, a = 0 and w = 6 are plotted as a function 
of Xr and \i. Here the constants used are II = 1 km, 
24 = 15 km, 2 c = 1 km and 26 = 0.4 km. The intersections 
of the zero Re(Sk) and zero Im(T^) lines give the 
eigenvalue solution X. It is seen that there exist many 
solutions with the most (and only) unstable solution 
having the smallest Xr or the longest wavelength. The 
solutions are labeled by an integer k with increasing k 
indicating increased \r or decreased vertical wavelength. 

The four longest vertical wavelength solutions are 
shown in Fig. 4 as a function of heating strength tn. 
It is seen that |Xil<<CXr and that the growth rates as 
represented by — X^ increase with m while the vertical 
wavelengths are almost constant. This diagram also 
indicates that the k=\ solution is most unstable 
throughout the range m = 4 to w=10, which cor- 
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Fig. -I. The four longest vertical wavelength solutions as a 
function of the heating rate m. The solid curves are Xr and the 
dash-dotted curves X^-. The corresponding vertical wavelength 
Z,t = 27r/Xr is given in the parentheses of the Xr coordinate. The 
ratio of static stabilities between stratosphere and troposphere 
is 1. 

responds to maximum heating rates of 4.7 and 11.7^C 
day”^ if u>6=0,29 cm s“h In fact, it is the only un- 
stable solution up to w = 9. For this solution the value 
of X for w = 5.8 is 0.388— 0.026i [km“^]], which is well 
within the range predicted by (21) and (23). 

The sensitivity of the eigenvalues with respect to 
the vertical heating profile (13) is shown in Fig, 5 in 
which the ^ = 1 to 1 = 4 solutions with w = 6 are given 
as a function of a. It is clear that the vertical wave- 
lengths are all insensitive to a. For the growth rates 
the k — \ solution is almost independent of a, but the 
other solutions all increase with a. There is an indica- 
tion that as a increased to >5 the shortest vertical 
wavelength becomes the most unstable solution. How- 
ever, for the realistic range of 1, which 

represents the statistical mean state of the heating 
function in the tropical atmosphere, 1=1 is the only 
unstable solution. 

Solutions for a different static stability in the 
stratosphere are obtained by including the radiation 
condition for z>Zt as the third boundary condition. 
In this case the reflection coefficient is derived by 
continuity requirements of w and <j> across Zt\ 

\-{Ss/StY^ 

r= exp(— 2 iXf 2 ;f), for Ss>St. 

1 + {S./Sd^ 

The result for Ss/St = 2,5 is shown in Fig. 6. It is 




Fig. 5. As in Fig. 4 except that the solutions are given as a function 
of the heating profile parameter a. 



seen that the general conclusions drawn from Fig. 5 
remain valid, although waves with a shorter vertical 
wavelength become more unstable than the k=\ waves 
beginning at a = 2 which corresponds to a ma.ximum 




Fig. 6. As in Fig. 4 except that the solutions are given as a 
function of the heating profile parameter a and that the ratio of 
static stabilities between stratosphere and troposphere is 2.5. 
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heating level of ^^12 km or 200 mb, much above all 
observed levels. So our result is basically insensitive 
to the vertical heating distribution. 

h. Vertical profiles 

The vertical structure of w, corresponding to the 
\ for k— 1 and a=0 in Fig. 6, is shown in Fig. 7 for two 
heating intensities ;w=6 and 10. It can be seen that 
the tropospheric amplitude reaches a maximum be- 
tween 6 and 9 km and has a minimum at the tropo- 
pause for both cases. In the stratosphere the amplitude 
for m = 6 increases at a slow rate with height due to 
the density factor, but the kinetic energy actually 
decreases with height. This decrease is due to the 
instability of the waves which results in Xi<0. The 
decrease of the wave kinetic energy is much faster 
for the more unstable w=10 case, so that even its 
amplitude decreases with height. The vertical phase 
distributions for both cases are almost the same and 
indicate that there is virtually no vertical propagation 
in the troposphere in contrast to a downward propa- 
gation in the stratosphere. The vertical wavelength 
computed from 27r/Xr is 16.2 km in the troposphere 
and 10.2 km in the stratosphere, because the value 




IiG. 7. Amplitude and phase of il; as a function of height for 
w = 6 (solid lines) and ;«=10 (dotted lines). The phase lines for 
both cases are approximately the same. 




of Xr is increased in the stratosphere by approximately 
a factor of as can be seen from (10). 

Vertical structures of other variables can be com- 
puted utilizing the basic equations (l)-(5) and the 
solutions to the meridional structure equation (6). 
Here only the temperature structure for the Kelvin 
wave case (;z= — l) will be shown. From (4) and (8) 
we have 




It is obvious from (25) that the amplitude of the 
temperature perturbation decreases as the zonal wave- 
number k increases. The Tn=~\ structure corresponding 
to the w in Fig. 7, is shown in Fig. 8. It can be seen 
that the amplitude distribution is quite sensitive to 
the strength of the heating. The discontinuity at Zt is 
due to the layered static stability distribution assumed 
in our model. The large temperature perturbations in 
the stratosphere are, of course, due to the large static 
stability there. Figs. 7 and 8 are constructed using 
nondimensional units. The two figures may be com- 
pared by letting the unit for the w amplitude be 
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! cm s“^; in this case the corresponding unit for the 
" amplitude will be a~^ where a = kX (radius of 

he earth) is the number of waves around the equator, 
"or these units the amplitude of w is quite close to 
i he typically observed values of “matured^’ tropical 
vaves, with maximum ^3 cm s~^ for m = 6 and 
^4.2 cm s~^ for m=10. On the other hand, the maxi- 
num T perturbation in the troposphere is almost 
10°C for wavenumber (a) 1 and even larger in the 
stratosphere. Lindzen (personal communication) has 
Suggested that these enormously large temperature 
.fluctuations indicate that large damping effect of 
(vertical momentum transport due to cumulus con- 
[vection must be important. However, for shorter waves 
I'the temperature perturbation is proportionally smaller, 
'if we consider a zonal wavelength of 4000 km, then 
^^^==10 and the temperature fluctuations are all <1°C 
in the troposphere, in good agreement with observa- 
jtions and with Holton (I971)’s numerical calculations. 
‘In fact, this magnitude estimate may explain why 
atmospheric Kelvin waves, which have a T amplitude 
in the lower stratosphere, are not observed in 
the troposphere. For such waves the vertical velocity 
in the troposphere would be an order of magnitude 
smaller than for the easterly waves which have a 
typical zonal wavelength ^3000 km and a maximum 
vertical velocity ^3 cm s“h On the other hand, damp- 
ing effect of a certain magnitude must play a role for 
the easterly waves to explain the complete absence 
of them in the stratosphere. Although these estimates 
are based on Kelvin waves onl\’, for other wave modes 
a relationship between the tropospheric temperature 
amplitude and the zonal scale must also exist at 
latitudes somewhat away from the equator because 
of the following scaling arguments. Wallace (1971) 
and Holton (1972b) have shown that, for weak, 
tropical, synoptic-scale motions with a depth scale 
comparable to the scale height, a diabatic heating 
rate of 5-lO^C day“^ is almost completely balanced 
by the adiabatic cooling. Thus the temperature fluc- 
tuation is one order of magnitude smaller than either 
term. This balance may be attributed to the smallness 
of the Rossby number which is '^O(l) or larger for 
tropical synoptic-scale motions. For planetary-scale 
motions the length scale is one order larger so that 
the Rossby number is '^0(10“^) at intermediate 
tropical latitudes, and the scale analysis indicates that 
the pressure and temperature perturbations must be 
increased by one order of magnitude. 

The phase diagrams in Figs. 7 and 8 indicate that 
in the lower troposphere and the entire stratosphere 
the temperature leads the vertical velocity by cycle. 
In the middle and upper troposphere (7-12 km for 
w = 6 and 5-12 km for w=10), they are in phase 
except near the tropopause (12-14 km), where they 
are out of phase. These phase relationships are con- 
sistent with the energetics because the conditions for 



wave growth are 



/■ 

/■ 

J 0 






Q'Vdz>^. 



From (25) and (10) it can be shown that these con- 
ditions are satisfied when Xi<0 and 



! SR^{iv)dz<Rl 
Jo Jo 



SRe(w)dz<R Qdz/iHcp), 

J 0 

The latter condition implies that the diabatic heating 
rate exceeds the adiabatic cooling rate which, of 
course, is responsible for the continuous growth of the 
waves with time. 



c. Equatorial trapping scales 

The meridional scale of the unstable waves may be 
estimated from the Hermite solutions of the meridional 
equation (6). The north-south, ^-folding width is 

>-e = ReC2(^/0V/3]s 

and the maximum amplitude for the mixed Rossby- 
gravity mode (/? = 0) occurs at 

yo=Re[(g/;)V/^]'- 

The maximum for the Kelvin mode occurs at the 
equator. From results obtained for the k = \ waves, 
\~0.388— 0.026z [km^^]] so that (^//)^'5=^ 25.5 — 1.7i (^m 
s~^] from (10). It follows that 1480 km and yo 
— 1050 km with /? = 2.21X 10~^^ m“^ s~K Both of these 
are reasonable values in terms of the trapping of 
waves on the equatorial beta-plane. 

6. Concluding remarks 

We have shown that tropical waves cannot remain 
unstable with respect to wave-CISK if the vertical 
wavelength is much smaller than the vertical scale 
of heating. This result holds for a number of reason- 
able heating profiles that have been observed or 
theoretically derived. The dispersive relationship of 
the growth rates of various equatorial beta-plane 
modes has been shown by Lindzen (1974) to be un- 
satisfactory in selecting the horizontal scale of tropical 
waves. From the computed eigenvalues it can be 
shown that, even under the extremely intense heating 
case of m=l0, the growth rates for Rossby waves, 
mixed Rossby-gravity waves, and the longer zonal- 
scale Kelvin waves are all very small. This result is 
similar to that obtained by Chang and Piwowar 
(1974) for the Ekman-CISK and indicates that wave- 
CISK is not responsible for the initial selection mecha- 
nism of the tropical waves. Once the waves are excited 
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by some other mechanism, however, the internal wave 
convergence could induce organized cumulus convec- 
tion and release wave-CISK as an energy source. The 
present calculation indicates that under such cases 
the shorter vertical-scale waves will still not be main- 
tained by wave-CISK. There is no definite explana- 
tion to the fact that the gravity modes and the shorter 
zonal scale Kelvin modes, which are usually not 
observed in the tropical atmosphere, are most unstable. 
However, these waves have very short periods and 
various arguments have been given against their oc- 
currence on a regular basis. For example, if one con- 
siders the lifting process^ of CISK (Ooyama, 1969; 
Holton, 1972b), the short periods will be insufficient 
for a low-level air parcel to reach its condensation 
level. In any case, the result that wave-CISK can 
neither excite nor maintain waves with short vertical 
wavelengths remains valid independent of these con- 
siderations. The concern about short vertical-scale 
excitation due to wave-CISK in a tropical numerical 
model may thus be relieved. The resonant forcing of 
tropical waves by a fundamental wave-CISK fre- 
quency associated with the zonally-symmetric, mixed 
Rossby-gravity mode as hypothesized by Lindzen 
(1974) also does not appear to be plausible, because 
the equivalent depth would be too large to give a 
fundamental periodicity near 5 days. 

The calculated vertical structure of w and T applies 
to all waves maintained by heating and is not re- 
stricted to wave-CISK. The eigenvalue X found by 
the wave-CISK equation used here may alternatively 
be viewed as being excited by a specified forcing 
[as in the model by Chang (1976)] with a small 
linear damping. The damping coefficient would then 
appear in the complex Doppler-shifted frequency in 
the same way an an imaginary frequency due to 
instability. The computed magnitudes of w and T 
suggest that the observed planetary-scale Kelvin 
waves are likely to be a result of forcing rather than 
instability, because a small forcing function can ex- 
plain the amplitudes in the stratosphere and tropo- 
sphere. However, an instability consistent with the 
present CISK parameterization should eventually 
amplify the waves to a finite amplitude w comparable 
to that observed for large-scale, cumulus-associated 
tropical motions, but with a temperature amplitude 
too large to be realistic. The small low-level vertical 
velocities of these waves is also an indication that 
CISK is unlikely to be important. On the other hand, 
the tropospheric amplitude of the synoptic-scale waves 
may quite possibly be supported by a CISK-type 
process. These waves have a wavelength 2000-4000 
km and a period ^5 days which lead to a phase speed 
relative to the ground of ~5-9 m s“L Such a phase 



2 Lindzen et al. (1975) have alternatively considered that the 
CISK is not a lifting process but a low-level organizing process 
which does not require lifting. 



speed would make the waves greatly attenuated before 
reaching the stratosphere, because the Doppler-shifted 
phase speed is very small near and below the tropo- 
pause and the damping can become very efficient 
(Lindzen, 1971). 
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ABSTRACT 

The problem of scale-selection of Kelvin waves in the stratosphere by forcing from tropospheric heating 
is analyzed using a simple linear model. The effect of vertical wind shear is excluded because the phase 
speed of the waves is fast relative to the range of the mean zonal wind in the vicinity of the tropopause at 
which level the upward energy flux due to forcing is evaluated. Results of this analysis modify Holton’s 
(1973) theory in that 1) the forcing is most efficient for the longest zonal wav^elength even if the heat sources 
are distributed randomly, and 2) the most favored vertical wav'elength of the excited waves is about twice 
the vertical scale of heating. The calculated vertical wavelengths e.xceed slightly those observed and the 
discrepancies are discussed. 



1. Introduction 

There is ample evidence on the existence of atmo- 
spheric Kelvin waves in the equatorial stratosphere 
(Wallace and Kousky, 1968; and others). These waves 
propagate eastward and downward with a zonal wave- 
number of 1-2, a periodicity of 10-15 days, and a 
vertical wavelength of 8-12 km. Their zonal wind 
component has a typical amplitude of 10 m s’”^ at the 
equator while no appreciable fluctuation of the me- 
ridional wind component is observed. The zonal wind 
and pressure perturbations are both symmetric with 
respect to the equator and are close to geostrophic 
balance in the meridional direction, but the waves 
behave like internal gravity waves in the zonal and 
vertical directions. More detailed observational and 
theoretical descriptions of these waves are given in 
Wallace (1973) and Holton (1975). These waves are 
believed to play a very important role in the quasi- 
biennial oscillation of the equatorial stratosphere by 
supplying westerly momentum to the zonal mean flow 
(Lindzen and Holton, 1968; Holton and Lindzen, 
1972). 

The vertical structure of the Kelvin waves indicates 
that they are maintained by energy propagating from 
the troposphere. In addition, the possibility that they 
are forced b\' interaction with mid-latitude motions, 
as suggested by Mak (1969), is quite small because 
of the small meridional wind perturbations which rule 
out the importance of lateral energy fluxes. Dynamic 
instability also appears unlikely due to the very long 
zonal scale and the fast westerly phase speed of the 
waves. 

Due to the presence of the large amount of deep 
cumulus convection in the tropical troposphere, it has 



been suggested that a plausible energy source fo 
these waves is the latent heat released by cumulu 
towers. There are two possible mechanisms by whic 
the latent heating in the troposphere may maintai 
the Kelvin waves: 

1) A two-way interaction such that the cumuli: 
convection, which supplies energy to the waves, i 
itself controlled, by the large-scale wave motion fielc 
This scale interaction results in an unstable situatio 
which is generally called the conditional instability c ■ 
the second kind (CISK). 

2) The waves are forced by cumulus heating whic 
is controlled by processes unrelated to the wave: 
Hayashi (1970) ’and Lindzen (1974) have investigate 
the first mechanism and found that the growth rate, 
of CISK are unsatisfactory in explaining the Kelvi 
waves. Holton (1972, 1973) studied the second mechr 
nism with a numerical diagnostic model in which th 
diabatic heating due to cumulus convection is spec 
fied. Damping in the form of Rayleigh friction an 
Newtonian cooling are also included. He found tha 
the wave response in the stratosphere to a trope 
spheric heat source resembling the Kelvin wave mod 
has a structure in close agreement with observations 
However, the observed periodicity of the Kelvi 
waves must be specified although no correspondin i 
peak in the cloud brightness spectra can be founii 
On the other hand, his numerical calculations sugge 
that the vertical scale of the excited waves is clos 
to that of the specified forcing. In equatorial wav, 
theory (Holton and Lindzen, 1968) the vertical seal! 
determines the wavenumber-frequency relationship 1 
By postulating a red-noise spectral distribution of th 
tropospheric heat sources, Holton suggested that th 
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vertical scale and the large low-frequency variance 
of the forcing result in a band-pass selectivity of the 
Kelvin waves. 

Holton^s theory may indeed be valid, especially 
because an oscillation of the monsoon heating in South 
Asia, with a quasi-periodicity of about 15 days, has 
been observed during the northern summer, although 
it is not clear whether this periodicity exists in other 
seasons. The purpose of this paper is to point out 
that the atmosphere acts like a band-pass filter even 
if the heat sources of the troposphere are distributed 
randomly in the frequency domain. The relationship 
between the vertical scale of the response and that 
of heating will also be examined. 

2. The forced model 

The zonal mean wind in the tropical stratosphere 
has considerable vertical shear which may influence 
the upward propagation of wave energy (Lindzen, 
1971). However, due to the fast phase speed of the 
Kelvin waves the shear is probably not very signifi- 
cant at the tropopause level for which our discussion 
will be most relevant. The zonal mean wind will 
therefore be neglected in our simple analysis. The 
effect of damping will also be excluded. The linearized 
zonal momentum, meridional momentum, hydrostatic, 
thermodynamic energy and continuity equations on 
an equatorial beta-plane, with a single zonal wave- 
number k, a Doppler-shifted phase speed c, and no 
perturbation in the meridional velocity may then be 
written 



--ikcu= -‘ikcpy 


( 1 ) 


d<t> 

l3yu= , 


( 2 ) 


dy 




d<f> RT 




dz h' 


( 3 ) 


Q 

-ikcT+wT = —, 


( 4 ) 


Cp 




d 






( 5 ) 



dz 

where w, T, <p and Q are the perturbation zonal 
velocity, vertical velocity, temperature, geopotential 
and the diabatic heating rate, respectively; H is o. 
constant scale height, T the static stability, Cp the 
specific heat at constant pressure, R the gas constant, 
13 the meridional gradient of the vertical component 
of earth’s vorticity, y the meridional coordinate, and 
z= -‘H \n{p ' po) is the vertical coordinate with p the 
pressure and pQ a reference pressure. 

Eqs, (l)-(2) give the meridional structure of the 



Kelvin waves: 



/ 

M, <^ocexp^-— j, (6) 



where the condition c>0 is required to satisfy the 
trapping condition on an equatorial /?-plane. If a 
height-dependence factor, exp( 2 / 2 //), is separated from 
the perturbation quantities, Eqs. (1) and (3)-(5) may 
be combined to form a single equation in w: 



where 



d'W 



4 


■XV = 


5 

2 


1 




AH- 


’’ = we~ 


•zKiH) 


R 

i' = 





z-2 ’ 



CpH 




( 7 ) 



( 8 ) 



The parameter X is a measure of the vertical wave- 
number. If the heating function is given, (7) can 
be solved with suitable boundary conditions. In our 
problem the following boundary conditions are used : 



jO, at 2 = 0 (9a) 

1 at z = Zt (9b) 

where 

Ci = rC2. 



Here Zt is the height of tropopause and the condi- 
tion (9b) results from the requirement that latent 
heating vanishes at Zi. The parameter r is a reflection 
coefficient which, in the absence of vertical wind shear, 
is given by the model specification of the static sta- 
bility distribution. The solution to (7) at z = Zt may 
now be written as 



w{zi) = / sin(\ 2 ) — dz^ (10) 

X(l+r) J zc 

where Zc is the height of the cloud base. 

3. The energy flux at tropopause 

We will consider the vertical energy flux 
where the overbar denotes the zonal ayerage and 
(^' = 0 exp( — 2 / 2/i), at the tropopause level ( 2 ,), to be 
a measure of the efficiency of forcing. Despite the 
actual existence of the vertical shear of the zonal 
mean flow in the stratosphere, the response of the 
waves should be largely determined by this quantity. 
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At Zt the solution is given by (9b) and only the Co 
term needs to be considered for upward energy propa- 
gation or downward phase propagation. 

Eqs. (1) and (5) may now be rewritten as 

— ikcu'= —{k(f>\ 

iku' — ( i\-\ Yy = 0, 

V 2hJ 

where u' = u exp(—z/2H). From these relationships it 
follows immediately that 




Since H^7 km for the atmosphere, \»l/(2//) if the 
vertical wavelength «88 km. This is the case for ail 
practical purposes so we may write 



cX 

4>' ^ — w' . 
k 



From the delinition of X given in Eq. (8) we may also 
write 



which leads to 



and 



cX«5q 

k 



— w'-. 
k 



( 11 ) 



Eq. (11) indicates that, if the amplitude of the ver- 
tical velocity perturbation at the tropopause level is 
h.xed, the upward energy flux will be inversely propor- 
tional to the zonal wavenumber. Examination of (7) 
reveals that icj' depends only on X (c is uniquely de- 
termined by X) as long as Q is independent of k. Thus 
for given vertical scale the forcing at the tropopause 
level becomes maximum for the longest horizontal 
wavelength. This result obviously remains the same 
if one considers the upward momentum flux u'w\ 
because u' = 4>'/c. 

In order to examine the selection of vertical scales, 
we assume a heating function 



mSq(y) sinTrf ), 

Q' = \ \ Sz / ( 12 ) 

lo, z>Zt or s<Zc 

where rn specifies the strength of the heating and 
Az = Zt-'Zc is a measure of the vertical scale of the 
heating. The function q(y) which specifies the me- 
ridional distribution of heating will be assumed to 



take the form of (6) for the time being. The vertic i 
profile as given by (12) resembles closely that of tl ' 
typical observations (Chang, 1976). For the simple ^ 
case, we shall first assume that the static stabililji 
parameter S has a constant tropospheric value 

5 = 5, 

throughout the atmosphere. In this case r = 0 and on 
upward energy propagation at z = zt is allowed. Tl. 
solution (10) at z, is now found to be 

w^=w exp[ — fXz,— /JyV (2^:)], z = s,, (L 

where 



7httAz X“+i/(4Hy rX "| /XAz' 

1 * 0 = sin -{zi+Zc) cos( 

X x2-X2Az2 L2 J \ 2 - 



(1^ 



denotes the amplitude of the solution in this simp 
model. We note here that no actual singularity exis 
in (14). The vertical energy flux at the equator ar 
z = z, is thus 

5»__ 

4>w' = — 

2k 



From (14) if one assumes X»l/(4/7"), it can be show 
that this energy flux reaches maximum at 



X 




The tropospheric forcing is therefore most efficient a; 
a vertical wavelength 

2x 

=2Az, 

X 

or twice the vertical scale of the heating. This resul 
is somewhat similar to those of Green (1965) an( 
Lindzen (1966). A plot of the energy flux profile at 2 ) 
is given in Fig. 1 in which it is assumed that z, = 15 2 ' 
and H = 7Zg. If we measure the vertical scale of Hoi 
ton^s (1972, 1973) heating function the value Az car 
be found to be ^^11 km, which predicts a maximuir 
response of 22 km vertical wavelength in the tropo 
sphere. Since in reality the stratospheric static sta 
bility is about three times that of the troposphere/ 
using (8) we may estimate that the stratospherk) 
vertical wavelength corresponding to a 22 km tropo! 
spheric value would be 

/5a-^ 

Ls^L,(^jj =22X3-*l«12.8km, 

where the subscript j denotes stratosphere. If we 
assume Az?^ 14 km, then 16.3 km. 



May 1976 



CHIH-PEI CHANG 



743 



It may be of some interest to consider the effect 
of different static stability between stratosphere and 
troposphere on the solution (13), because the ^‘effi- 
ciency” of reflection of wave energy at z = Zt may be 
scale-dependent. In this case we have 

S = Ss, for z>Ztf 
S = Sty for z^Zt. 

Another boundary condition is now needed and the 
simplest form available is the radiation condition 
above 2/, i.e., 

w' = A exp[— fX,2“‘/3yV(2c)], for 2 > 2 /, (15) 

where X, is the vertical wavenumber in the strato- 
sphere. The continuity conditions of and 0' across 
Zt lead to 

l-(5./5,)^ 

r= exp(~ 2 tXi 2 i), for Ss>Stf 

l + {Ss/Sd^ 

where X« is the vertical wavenumber in the tropo- 
sphere. The solution (10) at Zt in the form of (15) is 

2X< exp[i(\s — \t)zt2w 
(X,+X0~(X,-X0 exp{-2i\tZt)' 

where w is given by (14) with all X replaced by X«. 
The vertical energy flux at the equator and Zt is now 

5^ 

<t>W^—AA^^—Ew\ 



£ = 2[ (l+Ss/St) + ( 1 -S./Sd cos2X,2 

is the effect of different static stabilities above and below 
Zt and the asterisk denotes complex conjugate. For 
Ss—3Sty E varies between a maximum of 1 when 
Xf=0, Tv/zty lizlzty ..., and a minimum of 0.33 when 
X^=7r/(22/), 37T/(22i), 57r/(22i), .... Assuming 2^=15 
km and A 2 = 14 km, the variation of E is plotted in 
the upper section of Fig. 1 and the modified energy 
flux profile is indicated by the dashed curve in the 
lower section. It is seen that the most efficient forcing 
is slightly larger than twice the scale of heating. Hence 
we may conclude that the variation in static stability 
between troposphere and stratosphere does not cause 
great effect in the vertical scale of the most favored 
response. 

4. Concluding remarks 

We have shown that for randomly distributed 
tropospheric heat sources the Kelvin waves of the 





Fig. 1. Profile of vertical wave energy flux as a function of 
vertical wavelength Lz. The solid curve in the lower diagram is 
for the case S^ — St and the dashed curve for the case Ss = 3St- The 
upper diagram is the effect on E of the difference in static stability 
on the vertical wave energy flux profile for the case = 

largest zonal scale are excited most efficiently in the 
absence of damping. This may explain why most of 
the Kelvin waves observed in the lower stratosphere 
are of zonal wavenumber 1 and some are of wave- 
number 2. The vertical wind shear has been excluded 
from our analysis but it should not significantly alter 
this conclusion. 

The finding that the most favored vertical wave- 
length is about twice the vertical scale of heating 
gives a reasonable estimate for Kelvin waves in the 
lower stratosphere, although the calculated wave- 
length of 13-16 km is slightly larger than those ob- 
served. Several influences which are excluded in our 
model may alter this aspect. For e.xample, Lindzen 
(1971) found that the focusing effect of the vertical 
wind shear in the stratosphere would reduce the 
vertical scale of the waves. A two-scale expansion 
technique can be similarly applied to our forced model 
and the resultant prevailing vertical wavelength in 
the stratosphere should be reduced. 

Another effect which has been excluded is the 
actual north-south distribution of the heating func- 
tion. The meridional structure equation (6) indicates 
that the meridional scale of the waves is coupled wdth 
the vertical scale through the phase speed c. The 
variation of the meridional scale with the vertical 
scale may be considered slow because the e-folding 
width of the Gaussian distribution is proportional only 
to the square root of e, which is approximately pro- 
portional to the vertical wavelength. In the real 
atmosphere the meridional heating profile q(y) in (12) 
does not adjust to different vertical scales as we have 
assumed in our calculations. A proper treatment of 
this problem is, perhaps, to multiply the previously 
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calculated response by the projection 



f 

P= — 



q{y) exp[-^// (2c)]</j 



f exp[— /3// (2c)ydy 
J 0 
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Here we could specify q{y) as we did for the vertical 
heating profile, although there appears to be more 
observational uncertainty. For the purpose of the 
present qualitative discussion, however, this is not 
strictly necessary. We know that the latent heating 
is largely confined to the tropics, so that the meridional 
extent of q{y) must be quite limited. Therefore, the 
effect of (16) should be such that waves whose vertical 
scales are larger than that corresponding to the me- 
ridional scale of q{y) will suffer from having a large 
portion of their north-south domain unsupported by 
heating. If w^e assume that the meridional scale of q{y) 
can be represented by an e-folding width yc/~1500 km 
of a Gaussian distribution about the equator, and 
that 5e= 1.22X10”"* s”^ which corresponds to T = 3 K 
km“^, then those waves with a tropospheric vertical 
wavelength 

14.2 km, 

or a stratospheric wavelength 

Ls'^ 8.2 km, 

will have their amplitude reduced by (16). The reduc- 
tion will be more severe for longer vertical wave- 
lengths. This obviously would also shorten the pre- 
viously calculated most favored vertical scales of the 
Kelvin waves. 

Acknowledgments. The author wishes to thank Profs. 
J. R. Holton and R. S. Lindzen for discussion and 



REFERENCES 

Chang, C.-P., 1976: Vertical structure of tropical waves main- 
tained by internally induced cumulus heating. J. Atmos- 
Set., 33, 729-739. 

Green, J. S. A., 1965: Atmospheric tidal oscillations: An analysis 
of the mechanics. Proc. Roy. Soc. London, A288, 564-574 

Hayashi, Y., 1970: A theory of large-scale equatorial waves 
generated by condensation heat and accelerating the zonal 
wind. J. Meteor. Soc. Japan, 48, 140-160. 

Holton, J. R., 1972: Weaves in the equatorial stratosphere gener- 
ated by tropospheric heat sources. J. Atmos. Sci., 29, 368-375. 

, 1973: On the frequency distribution of atmospheric Kelvin 

waves. J. Atmos. Sci., 30, 499-501. 

, 1975: The dynamic meteorology of the stratosphere and 

mesosphere. Meteor. Monogr., 15, No. 37, 240 pp. 

, and R. S. Lindzen, 1968: A note on Kelvin waves in the 

atmosphere. Mon. Wea. Rev., 95, 385-386. 

, and , 1972: An updated theory for the quasi-biennial 

cycle of the tropical stratosphere. J. Atmos. Sci., 29, 
1076-1080. 

Lindzen, R. S., 1966: On the relation of wave behavior to source 
strength and distribution in a propagating medium. J. 
Atmos. Sci., 23 , 630-632. 

, 1971; Equatorial planetary waves in shear: Part I. J. 

Atmos. Sci., 28, 609-622. 

, 1974 : Wave-CISK in the tropics. J. Atmos. Sci., 31, 156-179. 

, and J. R. Holton, 1968: A theory of the quasi-biennial 

oscillation. J. Atmos. Sci., 25, 1095-1107. 

Mak, M.-K., 1969: Lateral driven stochastic motions in the trop- 
ics. J. Atmos. Sci., 26, 41-64. 

Wallace, J. M., 1973: General circulation of the tropical lower 
troposphere. Rev. Geophys. Space Phys., 11, 191-222. 

, and V. E. Kousky, 1968: Observational evidence of Kelvin 

waves in the tropical stratosphere, J. Atmos. Sci., 25, 900-907. 



.30 



Reprinted from Journal of the Atmospheric Sciences, Vol. 33, No. 8, August 1976 
American Meteorological Society 
Printed in U. S. A. 

Comments on “Instability Theory of Large-Scale Disturbances in the Tropics” 

C.-P. Chang 

Department of Meteorology^ Naval Postgrad^iate School^ Monterey, Calif. 93940 
4 February 1976 and 25 March 1976 



In a recent paper by Kuo (1975), a condition was waves and the Kelvin waves in the tropics. This condi- 
used in a CISK parameterization to explain the time tion, as stated by Kuo on p. 2231, is that 'The depletion 
and length scales of the observed mixed Rossby-gravity of moisture through precipitation must be replenished 

3 / 



1668 



JOURNAL OF THE ATMOSPHERIC SCIENCES 



Volume 3 



by evaporation from the surface during a complete 
period of the disturbance.” Eq. (15b) in his paper 
specifies this condition as follows : 

r (E+rate of moisture convergence) ^ maximum total 
precipitation given by deep cumulus convection 
during a period, 

where r is the period of the wave disturbances and E 
the evaporation rate. On p. 2234 he derived the limiting 
period ro for which the equal sign in (15b) holds. This 
is given by his Eq. (23) : 

Pi 

T0=— , 

E 

where Pi is the average maximum total precipitation 
over a period. This maximum total precipitation, 
according to the argument on p. 2232, is always pro- 
duced by the disturbance ‘‘once the convective insta- 
bility is activated by the large-scale flow,” because the 
“convection will proceed as an overturning process and 
produce the maximum amount of precipitation possible 
for the given air mass.” The value of Pt was estimated 
by subtracting the vapor contents of the outflow column 
from the inflow column as was schematically illustrated 
in his Fig. 3. The computed ro was then used as a cutoff 
period in the CISK parameterization [Eq. (24)] in such 
a way that the heating function is proportional to the 
factor r/ro it r<ro. The implication of his analysis is 
that, “disturbances with periods <ro cannot be sup- 
ported fully by deep cumulus convection initiated by 
such disturbances because the moisture content of the 
atmosphere will be depleted by the convective process.” 
The above arguments are, however, inconsistent with 
an essential part of the CISK parameterization used in 
the paper, namely that the large-scale heating rate is pro- 
portional to the vertical velocity W 5 at a low reference 
level Zh when IE6>0. In this parameterization when the 
wave amplitude at an early stage of the development is 
small, Wh is small and so is the heating rate, which 
implies that the precipitation rate is small. As the 
waves grow, Wh increases and the precipitation rate in 
turn increases. Since the total precipitation of a 



continuous wave chain cannot exceed the total evapora 
tion in the tropical belt, the restraining effect is sue 
that the precipitation rate must not e.xceed the evapora 
tion rate, or 3.75 g cm “2 day~^ according to the calcula 
tions quoted by Kuo. Therefore the total moistur 
availability in the tropics should constrain the maxi 
mum amplitude the waves can grow to, rather than th 
period of the waves. The intensity of the cumulns-scal 
convection is determined by the large-scale moistur 
convergence below in the CISK theory because thi 
convergence provides the high Se for the ascending air 
For more sophisticated parameterization theories, sucl 
as that of Arakawa and Schubert (1974) which is { 
form of the convective adjustment, the same incon 
sistency would occur. In these theories the effect 0 
cumulus convection is also entirely a response to tht 
large-scale forcing and therefore its intensity, as re 
fleeted by the net large-scale heating, should be con 
trolled by the large-scale flow. 

On the other hand, another condition mentioned b} 
Kuo [Eq. (15a)3 is a valid one, and may be quite 
important in selecting the wave scales. This conditior 
requires “a disturbance of sufficiently long duration and 
high intensity able to lift the lower atmosphere above 
the condensation level, thereby activating the deep 
cumulus convection.” In other words, waves with short 
periods are unlikely to grow. Kuo^s Eq. (24), which 
contains a cutoff period may, in a way, be viewed as a 
crude representation of this condition although the 
precise requirement can only be treated as a nonlinear 
problem. Based on this condition alone, waves whose 
Doppler-shifted phase speed equals zero for are 
most favored for development. This seems to be the 
case in the tropical western Pacific where the usually 
observed warm-core synoptic-scale waves have a 
typical easterly phase speed of 5-8 m s”*, quite close to 
the low-level mean zonal wind. 
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ABSTRACT 

In this paper we deal with the interpretation of observed oscillations in the tropical troposphere and 
stratosphere within the framework of the equatorial wave theory. A difficulty with this problem arises 
when one compares the short vertical wavelength (or equivalent depth) predicted by the classical theory 
and the observed large vertical scales associated with the low Doppler-shifted frequencies of the tropospheric 
oscillations. In this analysis it is shown that the inclusion of simple linear damping, justified by budget studies 
which revealed the important role of cumulus momentum transport, has a strong influence at low fre- 
quencies on the forced equatorial waves and results in two types of dispersive relationships. The first type 
is characteristic of the regular internal gravity waves which have fast phase speeds and weak vertical 
attenuation. The second type is dominated by the viscous damping time scale and has slow phase speeds 
and strong vertical trapping. The theory predicts that the stratospheric oscillations may be identified with 
the first type and the tropospheric oscillations with the second. In the case of Kelvin waves the results 
can be used to explain consistently both the observed stratospheric Kelvin waves and the planetary-scale 
Kelvin-like oscillations in the troposphere including the 40-50 day oscillation and the monsoon and Walker 
circulations. Possible implications with respect to other waves in the tropics are also discussed. 



1. Introduction 

Based on spectral and cross-spectral analyses of 
long time series (5^10 years) of tropical station data, 
Madden and Julian (1971, 1972) detected a global- 
scale oscillation in the tropical troposphere with a 
periodicity of 40-50 days. A degree of stationarity 
of this oscillation in time was also established by 
them by examining the limited amount of station 
pressure data available from the 1890^s. Their results 
indicate that this oscillation is most prominent in the 
zonal velocity component and the surface pressure, 
with the maximum amplitude situated along the 
equator. There are some time and spatial variations 
within a 40-50 day cycle as depicted in Fig. 1 which 
is reproduced from a schematic diagram constructed 
by Madden and Julian (1972), but basically the oscil- 
lation propagates eastward at a fairly steady phase 
speed and, at most times, exhibits a zonal wave- 
number 1 structure. Madden and Julian (1972) argued 
that enhanced deep cumulus convection is usually 
associated with large-scale upward motion and the 
oscillation appears as two vertical circulation cells 
along the equatorial zonal plane. 

The above structure resembles that of a theoretical 
atmospheric Kelvin wave (Holton and Lindzen_, 1968) 
in several ways : the direction of propagation, the 
meridional structure and the absence of large amplitude 
in the meridional wind component. However, for a 



given phase speed as slow as that observ'ed, Kelvin 
wave theory would predict a fairly short vertical 
wavelength which does not appear to agree with that 
of the observed oscillation. This discrepancy, which 
has been a major difficulty for the theoretical under- 
standing of the 40-50 day oscillations, has also been 
encountered by several meteorologists in their efforts 
to interpret various t>qDes of large-scale tropical wave 
motions using the equatorial wave theory. In a linear 
analysis Lindzen (1967) has shown that for synoptic 
(and larger) time scales, waves on an equatorial 
beta-plane generally propagate in the vertical as in- 
ternal gravity waves with rather short vertical wave- 
lengths. This result was evident in a study of the 
response of the tropical atmosphere to stationary 
heating by Webster (1972, 1973). In a two-level 
(750 and 250 mb) numerical model he found that the 
equatorial response in an easterly basic current has 
a horizontal structure that is characteristic of Kelvin 
waves. However, his analytical solution for the case 
of an idealized basic flow possesses very short vertical 
wavelength [~0(1 km)] and appears to suggest that 
higher vertical resolution is required for the nu- 
merical model. 

The short vertical wavelength predicted by the 
equatorial wave theory arises from the internal gravity 
wave character of these waves. Since for pure gravity 
waves the phase speed c is given by c={gh)^, where 
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Fig. 1. Schematic depiction of the time and space (zonal 
plane) variations of the circulation cells associated with the 
40-50 day oscillation. The mean pressure disturbance is plotted 
at the bottom of each chart with negative anomalies shaded. 
Regions of enhanced large-scale convection are indicated sche- 
matically by the cumulus and cumulonimbus clouds. The relative 
tropopause height is indicated at the top of each chart. (From 
Madden and Julian, 1972.) 

g is the gravitational constant and h a depth scale 
for the motion (usually called the equivalent depth), 
the vertical wavelength tends to be proportional to 



the phase speed. Within the range of large-scale 
tropical motions in the tropical troposphere the phase 
speed is usually slow compared to the external gravity 
wave phase speed, so the vertical wavelength is con- 
sequentially small compared to that given by the 
scale height. The observed Kelvin waves (Wallace 
and Kousky, 1968) and mixed Rossby-gravity waves 
(Yanai and Maruyama, 1966) in the lower strato- 
sphere have rather large Doppler-shifted phase speeds, 
(relative to the mean zonal wind) so that their ob-, 
served vertical wavelengths are in better agreement 
with the theory. In Table 1 we summarize the com- 
parison between the observed vertical wavelengths 
and that given by the equatorial wave theory for 
three major zonal wind oscillations in the tropics: 
the stratospheric Kelvin wave, the 40-50 day oscilla- 
tion and the planetary-scale stationary motion which 
could include the east-west overturning Walker cir- 
culation along the equatorial plane (Bjerknes, 1969; 
Webster, 1973) or the winter monsoon (Krishnamurti 
el al.j 1973; Webster, 1973). The theoretical vertical 
wavelength is calculated using the Doppler-shifted 
phase speed with the mean zonal flow specified by 
an averaged tropospheric easterly speed of — 5 m 
This specification is based on the evidence that all 
these oscillations, including the stratospheric Kelvin 
wave, are generated by energy sources in the 
troposphere. 

The equatorial wave theory used in Table 1 is 
basically inviscid because it does not include a sig- 
nificant frictional effect in the free atmosphere. In a 
numerical diagnostic study of the seasonal-mean cir- 
culation at 200 mb during the northern summer, 
Holton and Colton (1972) found that a very strong 
damping term (corresponding to linear damping time 
scale of ~1 day) is needed in the linearized barotropic 
vorticity equation to balance the generation of plane- 
tary-scale vorticity by horizontal divergence. Since 
the planetary-scale divergence field is largely a result 
of the strong condensation heating in the middle 
troposphere due to the vigorous cumulus convection 
associated with the summer monsoon, Holton and 
Colton hypothesized that the vertical mixing of 
vorticity and momentum by cumulus convection 
could account for this strong damping. This hypothesis 
is also consistent with several budget studies of syn- 
optic-scale, cumulus-related tropical motions (Reed 
and Recker, 1971; Williams and Gray, 1973; Reed 
and Johnson, 1974; Chu, 1976; and others), and the 
numerical diagnostic study by Chang et aL (1975). 
Recently Stevens et al. (1977) also show that this 
strong damping mechanism is responsible for reducing 
the theoretically calculated amplitude of the tropical 
waves forced by heating to the observed magnitude. 
Thus the inclusion of this mechanism in the equatorial 
wave theory seems to be necessary for the proper 
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Table 1. Comparison of observed tropical zonal oscillations and the in viscid equatorial wave theory. The theoretical vertical wavelengths 
are based on the observed Doppler-shifted phase speeds which assume a basic zonal wind of « = — 5 m s”^ See text for details. 









Observations 






Theory 


Oscillations 


Period 

(days) 


Zonal 

wavenumber 


Ground 
phase speed 
(m s-*) 


Doppler- 
shifted* 
phase speed 
(m s“‘) 


Vertical 

wavelength 

(km) 


Vertical 

wavelength 

(km) 


Stratospheric Kelvin waves 
(Wallace and Kousky, 


10-15 


1-2 


25-45 


30-50 


10-15 (strat) 
17-26(tropo) 


17-29 


1968) 

40-50 day oscillations 
(Madden and Julian, 


40-50 


1 


9-11 


14-16 


15-30 


8-9 


1972) 

Stationary planetary 
circulations (monsoon 
and Walker) 


QO 


1-2 


0 


5 


15-30 


3 



* M - 5 m s-^ 

interpretation of the observed tropical motions forced 
by cumulus heating. 

The purpose of this paper is to study the effect 
of the strong damping mechanism called for by pre- 
vious budget studies on tropical wave dynamics. The 
simplest form of linear drag will be used to represent 
this effect as was done by Holton and Colton (1972). 
We will show that this effect in some cases is more 
than just an expected reduction in the amplitude of 
the thermally forced waves or a near-trivial correction 
in the frequency equation. Moreover, in these cases 
it can account for most of the discrepancies in Table 1 
and provide a consistent theoretical basis for the 
interpretation of the tropical oscillations in both the 
stratosphere and the troposphere. 

2. Basic equations 

Since the vertical structure equation of all equatorial 
waves is identical to that for two-dimensional gravity 
waves, we will consider only the Kelvin wave case 
for the purpose of simplicity. Setting the meridional 
velocity zero we may write the linearized zonal mo- 
mentum, meridional momentum, hydrostatic, thermo- 
dynamic energy and continuity equations on an 
equatorial beta-plane, with a single zonal wave- 
number k and constant phase speec c relative to the 
basic zonal flow as 

— ikcn = -ik<j>--Du, ( 1 ) 

d<j) 

^yu= , ( 2 ) 

dy 

d<t> RT 
dz H ’ 



Q 

--ikcT+wT== DTy (4) 

d 

iku+e^^^ — =0, (5) 

dz 

where u, w, T, <t> and Q are the perturbation zonal 
velocity, vertical velocity, temperature, geopotential 
and the diabatic heating rate, respectively; .S' is a 
constant scale height, T the static stability, Cp the 
specific heat at constant pressure, R the gas constant, 
the meridional gradient of the vertical component 
of earth’s vorticity, y the meridional coordinate; and 
z= —H \n(p/pQ) is the vertical coordinate with p the 
pressure and a reference pressure. The damping 
effect is represented by the linear drag coeflScient D in 
the zonal momentum equation. A Newtonian cooling 
term with the same coefficient is included in the 
thermodynamic energy equation in order to facilitate 
the analysis. 

Eqs. (l)-(2) give the meridional structure of the 
Kelvin waves : 

u,<t> oc exp {—13^1 2c ) , (6) 

where 

c^c^iDf ky 

and the condition Re(c) = c>0 is required to satisfy 
the trapping condition on an equatorial /S-plane. If a 
height-dependence factor, exp(z/2H), is separated from 
the perturbation quantities, Eqs. (1) and (3)-(5) may 
be combined to form a single equation in w: 

bhjt/ 

uxV = — , 

az2 

3 b' 



( 3 ) 



( 7 ) 
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S 1 



w'=we^^^ (2^)^ 

R 

Q'^—Qe~zn2H)^ 



5= 



R 

— r 

H 



( 8 ) 



The parameter X is a measure of the vertical wave- 
number. If the heating function Q* is given, (7) can 
be solved with suitable boundary conditions. In our 
problem the following boundary conditions are used: 



0, at 2 = 0 (9a) 

2£;= - 

at 2 = zt (9b) 

where 

Ci=rC2. 



Here Zt is the height of tropopause and the condi- 
tion (9b) results from the requirement that latent 
heating vanishes at Zt. The parameter r is a reflection 
coefficient which, in the absence of vertical wind 
shear, is given by the model specification of the static 
stability distribution. The solution to (7) may be 
written in the form 



w{t) = / sinXz — dZy z^Zty (10a) 

X(l+r) J zc gh 

sinXz Q' 

w{z) = / — dz 

X(l+r)j^ gh X(l+r) 



X 



r . 

/ sinXz — dZy 

J zc ?,h 



Zt>z>Zcy (10b) 



sinXz Q' 

w{z)= / (r^*^^+e"*^^) — dZy Zc^z^O. (10c) 

\{l+r) J gh 



Here Zt is the height of cloud top and Zc the height 
of cloud base. Chang (1976) has used the above 
equations to show that, in the case of Z)=0, the 
vertical energy flux of Kelvin waves due to pressure 
work is 

w' y ( 11 ) 



and Lindzen, 1966), Chang (1976) was able to she 
that the Kelvin wave response given by the in vise 
theory agrees with the observed stratospheric Kelv 
waves. 

If we consider downward phase propagation (u 
ward energy propagation) only, Re(X)>0. For 
~7 km and vertical wavelength «88 km, Eq. ( 
may be approximated by 

<ax=6'i. (i:. 

In the inviscid case, D=0 and c=^c so that X is pu 
real. In the case of c=Cr+iCi where Cr=c> 

and Ci=D/k>^’y thus X becomes complex, i.e., X=» 
+fXi. Here X,* describes the exponential increa. 
(Xi>0) or decrease (X,<0) of the amplitude wit 
height. Eq. (12) now becomes quadratic and has fl 
following roots : 

S^±{S-^\rWy^ 

Cri.ii= >0, (!:• 

2\r 

-55±(5-4X,V)^ 

K.ii < 0 - 0 ^^ 

2ci 



The subscripts I and II are being used to denote th 
plus and minus roots, respectively. It is obvious tha 
Cri>Crij and | X^i | < | X,u | , hence the first solutio; 
(which shall be called Mode I) has a fast horizonta 
phase speed and can propagate more efficiently int 
higher levels while the second solution (Mode II) i 
slower and attenuates faster above the source region 
In addition, Eqs. (13) and (14) indicate that the wave: 
are excitable only if 




so that the effect of damping causes a short-wave 
cutoff in the vertical or a long-wave cutoff in the 
horizontal. 

If we take the inviscid limit of Z)— >0, Eqs. (13; 
and (14) reduce to 



^:ri = ^VXr 
X,i=0 
^ril =0 



(16)| 

(17)1 



where the overbar denotes the zonal average and 
(^>' = 0 exp( — z/2Z7). Eqs. (11) and (7) indicate that 
for this case the upward forcing of Kelvin waves is 
most efficient for long waves. Augmented by the fact 
that weaves with a vertical scale equal to twice that 
of forcing are most efficiently excited (Green, 1965 



Comparing (16) and (12) it is seen that Mode I cor- 
responds to the regular inviscid solution. Substituting! 
(17) into (l)-(5) leads to a trivial solution of vanisb'^ 
ing wave amplitude. Thus Mode I may be called the_ 
“regular mode’^ and Mode II the “viscous mode, 
because the latter is non-trivial only when 
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Fig. 2. Real Doppler-shifted phase speed Cr as a function of vertical wavelength Lt and damping 
time. For the convenience of application to other types of equatorial waves, the abscissa can also be 
designated by (g/Or*, where h is the equivalent depth. 



Fig. 2 shows the real phase speed Cr as a function 
of vertical wavelength Lz{ = 2Tr \r) and the viscous 
damping time scale \/D. For Z) = 0 (damping time 
scale the dispersion relationship of Mode I is a 
straight line of Cr^ Lz, typical of the internal gravity 
waves. In addition, the vertical coordinate (cr=0) 
represents the trivial solution of Mode II. For 
the two modes are represented by the quadratic 
curves. The effect of damping on Mode I is a rela- 
tively small modification with Lz increasing as Cr in- 
creases. The modification for Mode II is more sub- 
stantial and Lz generally decreases as Cf Increases. 
The two modes are separ^ited by the line Cr=Ci. For 
Mode I Cr>Ciy i.e., the inertial time scale is faster 
than the damping time scale. The reverse is true for 
Mode II. 

The imaginary vertical wavenumber X, is shown in 
I Fig. 3 as a function of vertical wavelength and damping 
time scale. Here the magnitude of X» decreases with 
' increasing Lz for Mode I and increases with Lz for 
j Mode II. The Cr — Ci curve again provides the parti- 
I tion of the two modes. 

From Figs. 2 and 3 it can be seen that, if a certain 
^ vertical wavelength is most efficiently excited by a 

S given forcing function, both modes will be excited 
for The regular mode will have fast zonal phase 

speed and less vertical attenuation, while the viscous 
, mode will have slow zonal phase speed and is more 
, trapped in the vertical away from the source region. 



3. Response of specified heating 

The Kelvin wave response of the tropical atmo- 
sphere can be calculated if the vertical heating func- 
tion in (7) is known. As in Chang (1976) we will 
use the vertical heating profile 

Q'cc5^(y) Sinsr^-— Zt'^z^Zc 

Q' = 0, z>Zt or s<Sc 

where the function j(y), which specified the meridio- 
nal distribution of heating, is assumed to take the 
form of (6). Solution (10) will be evaluated by speci- 
fying the vertical scale of the heating 2c 

= 14 km, Z) = 2.3X10“® s“^ (damping time = 5 days) 
and 55 = 35;, where 5s and St are the static stability 
for the stratosphere and troposphere, respectively. 

Since Mode I is less trapped in the vertical we ex- 
pect that it may be more prominent away from the 
source region. Using (10a) we display the vertical 
wave energy flux due to pressure work, at Zt as 
a function of vertical wavelength Lz and zonal wave- 
number s{=ka, where a is the radius of the earth), 
in Fig. 4. The cross-hatched area in the lower part 
of the diagram indicates the short vertical wavelength 
cutoff required by (15). The maximum response occurs 
at zonal wavenumber 1-2 and vertical wavelength 
30-35 kra. As expected, this result is similar to that 
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Fig. 3. Imaginary vertical wavenumber \ as a function of vertical wavelength and damping time. 



obtained by Chang (1976) except that the vertical 
wavelength is slightly longer and wavenumber 2 be- 
comes comparable in importance to wavenumber 1. 
This vertical wavelength range, according to Fig. 2, 
implies a Doppler-shifted phase speed ~50 m s"^, in 
reasonable agreement with the observ'ed stratospheric 
Kelvin waves. 

For the highly trapped Mode II the vertically- 
integrated tropospheric kinetic energy is considered 
an appropriate measurement for the response spectrum. 



Fig. 5 shows this quantity as a function of Lz and s. 
Here the zonal scale selection clearly favors the lowest 
zonal wavenumber, but all vertical wavelengths 
beyond the short-wave cutoff have about the same 
response. Thus the most prominent viscous mode ex-, 
cited should be of zonal wavenumber 1 and vertical 
wavelength >17 km. For vertical wavelength ~17 km 
the Doppler-shifted phase speed is 10 m s“^ accord- 
ing to Fig. 2. These values agree quite well with the 
40-50 day oscillation observed by Madden and Julian 



for MODE I 




Fig. 4. Vertical wave energy flux due to pressure work at the tropopause level, 
<t>'w'y as a function of vertical wavelength Lt and zonal wavenumber 5 for Mode I. 
Units are arbitrary. 
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K.E. Response for MODE H 




Fig. 5. Tropospheric kinetic energy response for Mode II as a function of vertical 
wavelength and zonal wavenumber. Units are arbitrary. 



(1972). However, if the forcing function has a white 
noise distribution in frequency the maximum Mode II 
response is still somewhat smaller than the maximum 
Mode I response. The absence of detectable 10-15 
day Kelvin waves in the troposphere therefore indi- 
cates that the thermal forcing is of a red noise nature. 
How’ever, this requirement of red noise distribution is 
not so much for the explanation of the stratospheric 
Kelvin waves, as suggested by Holton (1973); rather, 
it is for the explanation of the dominance of low- 
frequency oscillations in the troposphere. The 10-15 
day Kelvin waves in this case are of relatively higher 
frequencies and can become prominent only above the 
troposphere where the lower frequency viscous mode 
is greatly attenuated. 

The above results are summarized in Table 2 within 
the “non-resonant selection” category. The stationary 
waves related to the Walker circulation or the winter 
monsoon may, on the other hand, be viewed as due 
to resonant excitation by the stationary heat sources, 
because the forcing function has a pre*selected zero 
ground phase speed. By assuming a zonal mean wind 



of — 5 m s“^ in the source region the Doppler-shifted 
phase speed of this resonance is 5 m s“*h Referring 
back to Fig. 2, response is possible only in the Mode II 
regime with a vertical wavelength of ~27 km for the 
case of 5-day damping time. This vertical wavelength 
agrees much better with the observed monsoon or 
Walker circulation than Webster’s (1972) inviscid 
calculation of '^l km. In addition. Fig. 5 still indicates 
that zonal wavenumber 1 should be clearly preferred 
over shorter waves if the heat source has a white 
noise distribution in the zonal scales. In the real 
atmosphere the differential heating due to land-sea 
differences, as in the case of monsoon, or due to sea- 
surface temperature variations, as in the case of 
Walker circulation, is usually confined to lower zonal 
wavenumbers; so it is not surprising to observe the 
large zonal scale of these circulations. On the other 
hand, Egger (1976) reported that the response to the 
imposed heating anomalies in a linear model almost 
always favors the lowest wavenumbers regardless of 
the zonal scale of the anomaly imposed on the model. 
His experimental result is entirely understandable in 



Table 2. The most efficiently excited Kelvin waves predicted by the viscous theory. 







Vertical 

wavelength 

(km) 


Zonal 

wavenumber 


Doppler-shifted 
phase speed 
(m s"*9 


Ground 
phase speed 
(m s"‘) 




Mode I 


31-35 


1-2 


50-58 


45-53 


Resonant selection 




Mode II 
Mode I 


^17 

(not excitable) 


1 


15 


^10 


Non-resonant selection 




Mode II 


27 


1 


5 


0 



Damping time = 5 days. 

Vertical scale of heating = 14 km. 
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view of our result. This low-wavenumber preference damping term, i.e. 
for the thermally forced Kelvin wave response is 

related to (11) and the explanation has been given Q 

by Webster (1973) and Chang (1976). wV «£)r. 

Cj. 



4. Further comparison with the 40-50 day 
oscillation 

As mentioned previously, the two solutions of (12) 
are separated by the curve Cr=a. Since for Mode I 
the inertial time scale is faster than the damping time 
scale, it follows that the main balance in Eq. (1) is 
between the first two terms, i.e., 

—ikcu^ik4>. 

In this case u and <j> are in-phase and the wave is 
characteristic of the internal gravity waves. For 
Mode II the damping time scale is faster so that the 
main balance is between the last two terms, i.e., 

0= — ik<t>— Du. 

Thus u and 0 have a quarter-cycle phase difference 
which is characteristic of a highly viscous wave motion. 
If we examine the schematic diagram of the 40-50 day 
oscillation (Fig. 1) depicted by Madden and Julian, 
it can be seen that, for most of the categories, u and 0 
have a phase relationship between in-phase and 
quarter-cycle difference. The important effect of 
damping is thus clearly indicated. 

If the vertical scale is large, the main balance in 
the thermodynamic energy equation is between the 
adiabatic cooling and the diabatic heating. The residue 
is appro.ximately balanced by the advective term for 
Mode I, i.e., 

Q 

— ikcT^ wF 

Cp 



In this case T, w will be in-phase if diabatic heating 
exceeds adiabatic cooling and half-cycle out-of-phase 
if adiabatic cooling dominates over diabatic heating. 
In Madden and Julian’s (1972) observational study 
they found that the surface pressure is half-cycle 
out-of-phase with the 300 mb temperature and in-phase 
with the 100 mb temperature. Since vertical motion 
IS close to in-phase with surface pressure as shown 
by Fig. 1, it may be inferred that T and w are in-phase 
at 300 mb and half-cycle out-of-phase at 100 mb. The 
former level is in the upper troposphere where diabatic 
heating must dominate for the motion to be self- 
sustaining, while the latter level is near the tropopause 
where the diabatic heating vanishes but the vertical 
motion may overshoot through the cloud top. The 
thermal structure of the 40-50 day oscillation thus 
also appears to be consistent with the viscous theory. 

When damping is included, the meridional structure 
of the Kelvin waves as expressed by (6) contains a tilt 
of the wave axis because c is complex, i.e.. 




Here the first part of the right-hand side gives a 
Gaussian distribution in the meridional direction with 
a Gaussian folding with yc given by 



and r, w will have a quarter-cycle phase difference. 

For Mode II this residue is mainly balanced by the W / ' 
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Fig. 7. Comparison of the meridional phase tilt of the 40-50 day oscillation at 
140°E, 170°N, observed by Madden and Julian (1972), and the value of the most 
efficient Mode II response predicted by the viscous theory. 



This meridional trapping scale is identical for both 
the regular and the viscous modes for a given vertical 
wavelength. In the case of the most efficient Mode II 
response discussed in the previous section, \yo\ 
= 1500 km and agrees quite well with the observed 
40-50 day oscillation. The second part of the right- 
hand side of (18) gives rise to a periodic variation of 
the waves in y, and therefore a meridional tilt of the 
wave axis. It also implies a poleward phase propaga- 
tion for the Kelvin waves because Cr>0, Ci>0 and 

Xr>0, 

A strong meridional tilt was found by Madden and 
Julian for the 40-50 day oscillation using cross-spectral 
analysis of the surface pressure data. Their phase 
diagram, constructed by using the Canton Island series 
as the base series, is reproduced here in Fig. 6. Al- 
though the meridional phase tilt seems to vary in 
different geographical locations, near Canton Island 
and near the equator the tilt is primarily NW-SE 
north of the equator and NE-SW south of the equator. 
For eastward propagating waves this indicates a pole- 
ward phase propagation. In Fig. 7 the phase tilt of 
the 40-50 day oscillation deduced from Fig. 6 at two 
different longitudes, 170°W (near Canton Island) and 
140°E (active convection area in Indonesia), are 
plotted along with that given by the viscous theory. 
Within the limitation of our simple theor\^ and the 
approximations invoked, the theoretical result can be 
considered as fairly good in simulating the observed 
phase tilts. 



5. Concluding remarks 

We have shown that the inclusion of a damping 
term can result in two types of dispersive relationships 
for the forced Kelvin waves. The first type (Mode I) 
is characteristic of the regular internal gravity waves 
which propagate rapidly in the zonal direction and 
attenuate slowly in the vertical. The second type 
(Mode II) is dominated by the viscous damping time 
scale and has slower zonal phase speed and faster 
vertical attenuation. This viscous mode reduces to a 
trivial solution in the inviscid limit so its existence 
is supported only by damping. Assuming that the 
main thermal sources in the tropics are confined to 
the troposphere and have a red noise distribution 
which enhances the response of the viscous mode, 
this theory" predicts that the stratospheric waves 
behave like the regular mode and the tropospheric 
waves behave like the viscous mode. Our theory in- 
cludes features that are consistent with several results 
found by previous theoretical or numerical studies of 
the equatorial waves that incorporate the damping 
mechanism. Lindzen (1971) and Holton and Lindzen 
(1972) found that as a wave is propagating in weak 
vertical wind shear, the efficiency of damping would 
change inversely according to the change of the 
Doppler-shifted phase speed. This is consistent with 
our result that the slower moving waves are more 
trapped in the vertical. In a numerical diagnostic 
model of the Kelvin waves, Holton (1973) observed 
a maximum response in the tropospheric wavenum- 
ber 1 spectrum with period >30 days. This response 
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seems to correspond to our viscous mode. Very recently 
Stevens et aL (1977) reported that with a parameter- 
ization of the cumulus damping, the tropical waves 
have very little vertical structure. This is again ex- 
plainable by our theory because the vertical wave- 
length is increased as a result of damping. 

In real atmosphere the damping due to cumulus 
convection and other processes is much more com- 
plicated than the simple linear drag term used in our 
analysis. However, the agreement with observed zonal 
wind oscillations such as the stratospheric Kelvin 
waves, the 40-50 day oscillations, and the stationary 
monsoon and Walker circulations strongly indicate 
that our theoretical interpretation is quite relevant 
to these oscillations. In addition, although the present 
analysis is limited to Kelvin waves only, the vertical 
structure equation applies to all equatorial waves. 
Hence it would be quite possible that the general 
result also applies to other thermally controlled tropical 
wave disturbances. This is particularly interesting in 
view of the large vertical scale and slow Doppler- 
shifted phase speed usually observed in the tropical 
tropospheric waves such as the easterly waves. In a 
discussion of the present results, Wallace (personal 
communication) suggested that they w'ould also ex- 
plain the seemingly non-dispersive westward propaga- 
tion of the synoptic-scale cloud clusters observed in 
the time sections of tropical satellite pictures (Chang, 
1970). In these sections the propagating cloud patterns 
have quite uniform phase speeds but the apparent 
wavelength varies over a large range. This non-dis- 
persive character led Wallace (1972) to propose that, 
despite the strong vertical and longitudinal variations 
of the basic flow, they are more likely a result of 
advection by basic flow than a manifestation of some 
kind of Rossby waves. He argued that the difficulty 
of ascertaining a steering level is not as severe as that 
of explaining why the beta-effect does not produce 
dispersive phase speeds. However, to interpret these 
patterns as being passive cloud masses advected by 
the basic flow is also unsatisfactory due to the fact 
that they are frequently observed to be associated 
with active convections and that they sometimes 
decay and regenerate during the course of propaga- 
tion. These difficulties are partially alleviated in the 
framework of the present theory because the viscous 
mode should have a much smaller Doppler-shifted 
phase speed and would appear as if they are advected 
by the basic easterly current at some levels. 
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Some Theoretical Problems of the 
Planetary-Scale Monsoons 

By C tmi-Pia Cmano' ) 



Siimfuun Some importani theoretical problems of the planetary-scale monsoons which have arisen 
from recent advances of observational studies are reviewed. These include. ( 1 ) the requirement of a strong 
damping mechanism in the planetary scale vorticity budget of summer monsoon and a similar but weaker 
requirement for the winter monsoon: (2) the localized barotroptc instability of the summer monsoon 
w hich IS a result of the strong zonal asymmetry of the planetary-scale tlow and causes significant nonlinear 
energy conversions; and (3) the oscillations of the planetary-scale monsoons. It is pointed out that these 
problems are inter-related and their understanding is also important for the proper simulation of other 
scales of motion of the monsoon circulation. 

Key words: Monsoon; Oscillation, vorticity budget: Barotropic instability. 



/. Iniroduction 

The monsoon circulation of both the northern summer and the northern winter 
(hereafter the word 'northern' will be neglected) are characterized by large planetary 
scale quasi-stationar\ wa\es in the subtropics and tropics. As have been shown by 
Krishnamurfi (1971a), Krishnaml rti. Kanamitsi'. Kass and Lee (1973b) and 
others, these vva\es are a result of the thermally direct overturnings induced by the 
diiferemial heating associated with the land-sea distribution. The major component 
of the heating is the convective latent heat release which is most intense over the 
Tibetan and Mexican highlands during summer and the maritime continent of 
the Malaysia-Indonesia region during wanter. Mainly due to the availability of the 
commercial aircraft data, the principal features of these zonal asymmetries have been 
best depicted by the wind analyses near 200 mb by Krishna.ml rti and Rogers ( 1970), 
Krishnamurti et ai (1973b) and Sadthr (1975). During summer these features 
include a strong anticyclone centered over the Tibetan plateau, a weaker anticyclone 
over Mexico and elongated troughs oriented northeast-southwest in the mid-Pacific 
and mid-Atlantic. These are all well presented in Fig. 1 which shows the time-mean 
streamfunction at 200 mb for June August 1967 as determined by Krishnamurti 
( 1971a). During winter a large streamfunction high over the equatorial South China 
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Summer mean (June -August 1967) 200-mb streamfunction, contour interval is 50 ^ 10*^ m’ sec (f'rom 

Krishnamurti, 1971a.) 




Figure 2 

Winter mean (December February 1971) 200-mb streamfunction. contour interval is 50 • 10^ nr sec~F 

(From Krishnamurti er 1973b.) 



Sea and western Pacific replaces the Tibetan high to become the most pronounced 
planetary scale feature, with two weaker highs appearing over the equatorial eastern 
Africa and the northeastern South America and three oceanic troughs e.xtending 
from the equator to the southern subtropics. These features are shown in Fig. 2 which 
is the 200 mb time-mean streamfunction for the winter of 1971 determined by 
Krishnamurti et aL (1973b). 

A number of interesting theoretical problems have arisen from the observed 
monsoon circulations. In this paper we will review some of the problems that are 
related to the planetary-scale monsoon. They include: 

(1) The requirement of a strong damping mechanism in the planetary scale 
vorticity budget. 
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(2) The localized barolropic instability associated with the planclary-scalc 
summer monsoon, and 

(3) Oscillations of the planetary-scale monsoons. 



2. The requirement of a stronq planetary scale clamping mec ham.sni 

In a linear, one-level numerical model, Holton and Colton (1972) used the 
barotropic vorticity equation to diagnose the budget of the seasonal-mean planetary 
scale vorticity at 200 mb during summer. They used the 200 mb mean divergence 
observed by Krishnamurti (1971a, b) as the steady forcing function with a mean 
zonal flow also specified by the observed values. Figure 3 shows the velocity potential 




Summer mean (June -August 1967) 200-mb velocity potential, contour interval is 10 10^^ m’ sec V The 

few streamlines sketched indicate the direction of the mean divergence motions. (From Krishnamurti. 

1971a.) 



computed by Krishnamurti (1971a) which expresses the observed mean divergence 
field. Comparing this figure with Fig. 1, it is evident that the time-mean vorticity and 
divergence fields are approximately half-cycle out of phase. However, the resultant 
steady-state planetary scale vorticity field calculated by Holton and Colton is 
about one-quarter cycle out of phase from the divergence field, unless a very strong 
damping term is added to the mean vorticity equation. This difficulty arises because, 
in the absence of strong damping, the main balance in the vorticity equation is be- 
tween the generation by horizontal divergence and the advection by mean flow. This 
near-balance requires that the vorticity be one-quarter cycle out of phase from the 
divergence. The very efficient generation by the strong divergence also causes the 
amplitudes of the long waves to be significantly larger than those observed. Figure 4 
(from Holton and Colton, 1972) compares the diagnosed and observed amplitudes 



1092 



Chih-Pei Chang 



(Pageoph, 



and phases of the streamfunction for zonal wavenumber one as a function of latitude. 
When the linear damping coefficient is increased from a modestly strong value of 
5 X lO'^sec'^ to a very strong value of 1.5 x 10“^sec“\ the discrepancy dis- 
appears and the phase and amplitude of the long waves become in close agreement 
with observations. 

PHASE (longitude) 



90*W 0“ 90'E I80' 90*W 




AMPLITUDE 

Figure 4 

Amplitude in units of 10^ sec" ‘ (left-hand curves) and phase (right-hand curves) of the streamfunction 
for zonal wavenumber 1 as a function of latitude. Solid curve is observed streamfunction: short and long 
dashes are the solutions computed by Holton and Colton (1972) for damping rate of 1.5 10"^ sec" ‘ 

and 5 ^ 10 ^ sec” respectively. (From Holton and Colton, 1972.) 

Holton and Colton's (1972) findings suggest that a very strong sink with a 
damping time scale faster than 1 day is required to balance the linear vorticity budget 
of the seasonal-mean 200 mb planetary-scale fiow during summer. This damping is 
much stronger than the normal dissipation process in the upper troposphere whose 
damping time is on the order of weeks. In looking for an explanation of this excessive 
damping requirement, they suggested that two possible mechanisms not included in 
their model may be responsible for the strong damping. The first is the vertical mixing 
of vorticity by cumulonimbus convection which is very strong over the Asian conti- 
nent and Mexican highland during the northern summer. This mechanism has been 
suggested by Riehl and Pearce ( 1 968). Reed and Recker (1971), Wallace (1971), 
Williams and Gray (1973) and others to balance the vorticity budget of the synoptic- 
scale easterly waves in the tropical troposphere. However, Holton and Colton 
(1972) felt that it alone is still insufficient to account for the required fast damping. 
They then proposed the second mechanism which is related to the day-to-day 
fluctuations of the monsoon vorticity and divergence fields observed by Krishnamur- 
Ti ( 1971b). Considering deviations from the time mean, the complete divergence term 
in the time-mean vorticity equation is <(C /)^> = (<C> T where C 

and / are the relative and planetary vorticiiies, respectively: S is divergence and the 
primes are deviation from the time mean < >. The last term, which was excluded in 
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Holton and Colton's linear calculations, tends to be negative according to the 
observed out-of-phase correlation of the fluctuating vorticity and divergence fields. 
This term therefore may counteract the vorticity generation by the time-mean 
divergence. Holton and Colton then speculated that the combination of this 
fluctuation correlation and the cumulus mixing mechanisms could be responsible for 
the less than one day damping lime. 

The requirement of a strong vorticity sink has also been confirmed by other 
numerical investigations that do not include the strong damping. As a follow up of 
the work of Holton and Colton (1972), Colton ( 1973) used a semi-spectral repre- 
sentation of the one-level nonlinear barotropic vorticity equation to study the scale 
interactions of the 200 mb summer monsoon flow forced by the observed mean 
divergence. As a consequence of the inclusion of only modest damping (time scale 
5 days), the resultant planetary-scale waves have amplitudes considerably larger than 
the seasonal mean values. Westward phase shifts of 50° for the Tibetan high and 30° 
for the mid-Pacific trough, in comparison with the observed positions, are also 
produced by the model (Fig. 5). Abbott ( 1973) developed a 3-level nonlinear quasi- 
geostrophic hemispheric model to simulate the time-mean climatology of the summer 



Streamfunction of the forced motion averaged over days 15 through 40 of the numerical integration by 



monsoon. To drive the monsoon circulation, he specified a steady heating distribution 
which was found by inverting the governing equations with the vorticity field giving 
by the observed time-mean values. This heating distribution {Fig. 6) in turn success- 
fully simulates the principal features of the planetary-scale flow as well as the nonlinear 
energy e.xchanges. However, it is about a quarter cycle out of phase from that implied 
by the observed mean divergence field and the satellite cloud data. This apparently is 
also due to the absence of a strong damping mechanism in the model which forces a 
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Figure 5 



Colton ( 1973). Contour interval is 50 ^ 10^ m“ sec ^ 
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Figure 6 

Time-independent horizontal distribution of the diabatic heating rate at 350 mb as computed by Abbot 

(1973). Isopleths are at 5 x 10 ' m“ sec " ‘ . 

balance between the advection and generation of vorticily. In a 3-component semi- 
spectral primitive equation model, Murakami (1974) also used a specified healing 
function to study the planetary scale monsoon flow during summer. Without the 
inclusion of the strong damping and allowing the specified heating to fluctuate with 
time, he was able to obtain the proper amplitude and phase of the long waves. 
However, the resultant divergence field in his model is somewhat diflferent from that 
observed by Krishnamurti (1971a). The greatest discrepancy occurs at 200 mb, 
where Krishnamurti (1971a) found very strong divergence in the Tibetan high area 
during the 1967 summer, while in Murakami's ( 1974) model this level is almost non- 
divergent. The model also produces a different horizontal distribution of long wave 
divergence compared to the observed 1967 distribution, with the Tibetan divergence 
center shifted slightly to the east and the mid-Atlantic convergence center shifted to 
the longitude of central Africa. Thus the magnitude of divergence forcing is much 
reduced and a phase shift between divergence and anticyclonic vorticity also results'). 
In addition, the advective effect is not properly treated because the zonal mean flow 
is not allowed to interact with the three interacting long waves. Hence it is unclear 
what is the main mechanism that produces the correct amplitude and phase of the 
vorticity field. 

Looking for alternative sources for the strong damping, Barcilon and Cooper 
(1974) hypothesized that a thermal boundary layer effect could operate in the 
atmosphere during summer monsoon. In their simple analytical model a sinusoidal 
temperature variation is applied at the lower boundary which provides the planetary 
scale differential heating. A number of solutions were found for a balance between 



■) However. Ml rakami (personal communication) found that this horizontal distribution of diver- 
gence is in good agreement with the 1970 summer data. 
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the horizontal advection and vertical diffusion of heat in the thermodynamic energy 
equation, and one of them allows a half-cycle out of phase shift between the upper 
level vorticity and divergence. However, since the monsoonal differential heating in 
the real atmosphere is mainly due to the latent heat release of cumulus convection 
which is maximum in the middle-upper troposphere, it is uncertain whether Barcilon 
and Cooper's (1974) theory is completely applicable to the real atmosphere. 

Another explanation for the observed vorticity-divergence phase relationship was 
offered by Paegle and Paegle (1976). They speculated that a high pressure region 
with strong pressure gradients may cause the geostrophic vorticity to become less 
than one-half of the negative planetary vorticity. This implies that the balance 
equation becomes non-elliptic and, under certain rather strong assumptions, leads to a 
rapid vorticity adjustment with strong horizontal divergence occurring in the negative 
relative vorticity region. However, the resultant steady-state absolute vorticity 
vanishes and this efficiently eliminates the generation of vorticity by divergence. Such 
large negative values of relative vorticity are not observed in the large-scale monsoon 
How. In fact the minimum seasonal mean relative vorticity at 200 mb during the 1967 
summer as observed by Krishnamurti is ^ —3 x 10~ ^ sec” ^ which occurs in the 
vicinity of the Tibetan and Mexican highs near 30°N. This magnitude is considerably 
smaller than that of the Coriolis parameter. Thus it appears unlikely that this super- 
gradient effect is relevant to the planetary-scale monsoon problem. 

To test the hypotheses proposed by Holton and Colton (1972), Chang and 
Pentimonti ( 1977) used a 3-level primitive equation model with a 4° finite-difference 
resolution in both the zonal and meridional directions. The zonal mean flow is 
included usinga Newtonian term which maintains the proper north-south temperature 
gradient by continuously adjusting each temperature toward a thermal equilibrium 
determined by the zonally symmetric climatology. The zonally asymmetric motions 
are generated by a specilicd perturbation heating function whose horizontal distribu- 
tion follows that of the 200 mb seasonal mean divergence observed by Krishna.murti 
(1971b). It turns out that this heating distribution almost reproduces the observed 
mean divergence in the upper troposphere, which suggests that the choice of the 
forcing function is proper. Two experiments were carried out by integrating the model 
in time until 25 days after it reached the respective steady states. Both experiments 
include only weak damping (time scale >7 days) in the free atmosphere but the 
heating function was held steady only in the first experiment and allowed to fluctuate 
in time with an amplitude equal to one-half of the time-mean heating in the second 
experiment. The fluctuation period was set to 10 days. Comparing the results of the 
two experiments, Chang and Pentimonti found essentially no difference between 
the time-mean planetary scale fields except a slightly reduced amplitude for the 
fluctuating case. Figure 7 is their time-mean wind field at 250 mb (near the upper level 
of the model) for the fluctuating heating case. It can be seen that the simulated 
Tibetan and Mexican anticyclones are shifted about 30® and 25®. respectively, from 
the observed positions. These phase shifts are almost identical to those produced by 
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Time mean wind field at 250 mb forced by a fluctuating heating function as computed by Chang and 
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their steady-heating case. Thus the fluctuating component as specified in the model 
is almost inconsequential in providing the damping needed for the phase shift prob- 
lem. On the other hand, in both experiments the amplitudes of the long waves are not 
excessively large compared to observations. This result is similar to that of Abbott 
(1973) and seems to suggest that nonlinear processes could be responsible for part of 
the planetary'-scale vorticity sink deduced by Holt( )N and Colton's linear calculations. 
The fluctuation of the planetary-scale flow apparently enhances these nonlinear 
conversions resulting in the further reduction of its time-mean amplitude. 

A noteworthy feature in Chang and Pentimonti's ( 1977) time-mean fields is that 
the positions of the upper tropospheric mid-oceanic troughs do not have significant 
phase errors as those of the anticyclones. As pointed out by Krishnamurti, Daggu- 
PATY. Fein, Kanamitsu and Lee (1973a), the damping due to cumulus transports 
should not operate over the relatively cloud-free oceanic areas. With the inclusion of 
nonlinear processes (and also large-scale vertical advections), the damping require- 
ment implied by the phase errors of Chang and Pentimonti's time-mean upper level 
anticyclones is weaker than that found by Holton and Colton (1972). Perhaps a 
damping time longer than 1 day, which seems more reasonable than Holton and 
Colton's value if interpreted as due to the cumulus transports, may be sufhcienl. 
In this case a numerical experiment similar to Chang and Pentimonti's, but in- 
corporating a moderately strong damping term with time scale order of a few days in 
the continental convccliv e area only, may prox ide some insight as to whether cumulus 
transport processes are responsible for the missing vorticity sink. 

The studies mentioned above all concern the summer monsoon only. During 
winter Krishna.viorti e! al. (1973b) observed that the rotational component of the 
seasonal-mean flow of the actixe monsoon region is almost one order smaller than 
that of the summer monsoon, xvhile the divergent component has a comparable 
magnitude. The reason for this appears to be the proximity of the winter monsoon 
convections to the equator, xvhich reduces the etficiency of vorticity generation by 
horizontal divergence. Hoxvcver, the importance of a dissipation mechanism in the 
winter monsoon was also brought out in a fexx theoretical studies recently. In a linear 
study of the stationary forcing in the tropical atmosphere, Webster (1972. 1973) 
found that the primary response in the equatorial region corresponds to the atmos- 
pheric Kelvin waxes (Holton and Lindzen. 1968). A theoretical calculation shows 
that the zonal wavelength of this Kelvin xvave response is very long while its vertical 
xvavelength is very short [0(1 km)]. The former agrees quite xvell xvith observed 
equatorial stationary circulations such as the xvinter monsoon, but the latter is one 
order of magnitude shorter than the observed vertical wavelengths. This discrepancy 
arises because the equatorial Kclx in waxes behave like internal gravity xvaves in the 
zonal-xertical plane and their phase speed, c, is related to an equivalent depth scale h 
by the relationship c = {(/h)^ “. These stationary xvaves are imbedded in a mean 
easterly zonal current so that they have a Doppler-shifted phase speed of a few meters 
per second, xvhich implies a x-ery short vertical wavelength. To get around this problem 
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in simulating the observed stationary motions, Webster (1972, 1973) adopted a two 
level numerical model with friction which is quite successful in reproducing many of 
the observed features. 

Very recently Chang (1977) re-examined this problem by incorporating linear 
damping terms in the thermally forced Kelvin wave equations. The dispersive rela- 
tionship of the Kelvin waves may be written as ck — where k is the vertical 
wavenumber and S is a static stability parameter. The linear damping introduces an 
imaginary part to c which causes k to become complex also. Consequently the 
dispersive relationship becomes quadratic and has two types of solutions. The first 
type arises when the inertial time scale is faster than the damping time scale. This 
solution resembles the regular internal gravity waves whose vertical wavelength 
decreases as the Doppler-shifted phase speed decreases. The second type, called the 
viscous mode by Chang ( 1977), arises when the damping time scale dominates over 
the inertial time scale and its vertical wavelength increases as the Doppler phase speed 
decreases ( Fig. 8). This solution is also much more trapped in the vertical so that it has 
appreciable amplitude only in the region of thermal forcing (troposphere). In the 
inviscid limit this solution reduces to a trivial solution with vanishing amplitude. For a 
reasonable range of physical parameters, Chang (1977) calculated a vertical wave- 
length of > 17 kms for the stationary response, which is a significant improvement 
over the inviscid theory. The horizontal structure of this solution is also in good 




Real Doppler-shifted phase speed. as a function of vertical wavelength. L,, and damping tune for the 

Kelvin waves. (From Chang. 19"7.) 
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agreement with observations and Webster’s (1972, 1973) numerical model. This 
requirement of a damping mechanism with time scale on the order of 5 days is much 
weaker than that found by Holton and Colton (1972) for the summer monsoon 
and could well be explained by cumulus transports. Detailed observational studies of 
the winter monsoon are needed to improve the theoretical understanding of this 
problem. 



3. Localized hctrotropic instability of the planetary scale summer monsoon 

As mentioned previously, the amplitude of stationary long waves during summer 
is better simulated by the nonlinear models than by the linear models when excessively 
large damping is excluded. This apparently is due to the barotropic energy conversions 
from planetary to both zonal and synoptic scales which are allowed in the nonlinear 
models only. Observational evidence of this type of energy conversion has been 
presented by Kana.mit.su. Krishnamurti and Depradine (1972). who computed the 
wave- zonal flow and wave-wave energy exchanges for the tropical belt from 15°S to 
15°N using the 1967 summer 200 mb wind data analyzed by Krlshnamurti and 
Rogers ( 1 970). In particular, they showed that the synoptic scale waves receive kinetic 
energy at a fairly significant rate (order of lO^^ m- sec '^) from both the zonal flow 
and the planetary-scale waves. A similar result for the equator to 30°N belt was also 
obtained by Krish.namurti et al. (1973b). 

Transient synoptic-scale wave disturbances have been observed in the tropical 
upper troposphere in all seasons. During summer, the most active season, the distur- 
bances are found to occur primarily in two types of regions. One type is in the vicinity 
ofthe mid-oceanic troughs (Rick. 1959: Sadler. 1967: Miller and Carlson. 1970: 
Frank. 1970) where the disturbances may assume the form of open waves or closed 
vortices. Occasionally these disturbances may develop in the vertical or be coupled 
with low level waves, and intensification leading to tropical cyclone development 
sometimes may also follow. Another region of transient wave activity is alone the 
upper tropospheric easterly jet over the Indian ocean, where Krishna.murti and 
Rogers ( 1 970) have reported fairly regular occurrences of synoptic-scale disturbances. 
Apparently these regions are where most of the barotropic energy conversion from 
zonal and planetary scales to the synoptic scale takes place. Thus it appears that the 
nonlinear energy transfers of the summer monsoon have a persistently localized 
nature. Another evidence of this is related to the linear barotropic instability criterion, 
which has been shown by Merilees (1968) to be usually satisfied if nonlinear energy 
e.xchanges are very prominent. The upper tropospheric zonal mean flow of the 
summer monsoon does not satisfy this criterion, but regions of large negative 
gradients o( relative vorticity (Fig. 9) where this instability criterion is satisfied 
locally can be found in the seasonal mean 200 mb analysis (Krishna.murti. 1971b). 
These negative vorticity gradients occur mainly in the vicinity of the easterly jet and 
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Figure 9 

Summer mean (June-August 1967) relative vorticiiy at 200 mb, contour interval is 10'^ sec“‘. (From 

Krishnamurti, 1971b.) 



the mid-oceanic troughs, and on some individual days their magnitude may even 
greatly exceed the seasonal mean values. Therefore it is plausible that at least some of 
the transient synoptic wave activities in these regions is generated by the localized 
(and sometimes short term) barotropic instability, which is clearly a result of the 
zonal asymmetry of the planetary-scale summer monsoon. 

The kinetic energy conversions from the planetary scale and the zonal flow to the 
synoptic scale in the upper troposphere have been simulated numerically by Abbott 
(1973) and Krishnamurti ei al. (1973a). The conversions produced by Abbott’s 
(1973) model are roughly of the same order of magnitude as thos^ computed by 
ICanamitsu et uL ( 1 972). However, the conversions obtained by Krishnamurti et al, 
(1973a) are about one order smaller. This appears to be related to their treatment of 
the Tibetan high as a fixed barrier which prevents vorticity transport across the 
boundary of the barrier. As a consequence, barotropic instability cannot be realized 
in the immediate vicinity of the Tibetan high and the total amount of energy extracted 
from the larger scales by the synoptic scale within a fixed period of time must be 
greatly reduced. Nevertheless, Krishnamurti et aL (1973a) were able to simulate the 
gross features of the planetary scale waves in the upper troposphere, such as the 
mid-oceanic troughs and the Mexican high. Apparently once the strength, shape and 
position of the Tibetan high are fixed in the experiment, the interaction of the zonal 
mean flow and this blocking high is sufficient to produce mechanically many other 
features of the circulation. 

In the numerical model by Colton (1973), attention is focused on the development 
of the synoptic-scale disturbances in the easterly jet and in the mid-Pacific trough 
regions where the barotropic instability criterion is locally satisfied. Near the easterly 
jet the model produces disturbances that propagate westward at a phase speed slightly 
less than the local speed of the Jet and continue to grow until they are well past the 
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longitude of the jet maximum. After reaching the maximum amplitude they begin to 
decay fairly rapidly. Many features of these disturbances resemble those analyzed by 
Krishnamurti and Rogers (1970) in the Indian Ocean and African region. In the 
mid-Pacific trough the model also produces waves and vortices whose movement and 
behavior are quite similar to those observed in daily synoptic charts. In both regions 
the growing disturbances exhibit a latitudinal tilt in opposition to the lateral shear of 
the quasi-steady basic flow, which is characteristic of the barotropically unstable 
waves. 




Figure 10 

Zonal average kinetic energy at 20N as a function of time and wavenumber as computed by Colton ( 1973). 
Contour interval is 5 m' sec'C Hatched areas are regions of energy maximum. 



In studying the barotropic energy exchanges, Colton (1973) also examined the 
time evolution of the computed kinetic energy spectrum. An example of this is given 
in Fig. 10 which shows the zonally averaged energy at 20'^N as a function of time and 
wavenumber. A noteworthy feature is that when the energy of wavenumber 1 is at a 
maximum, the energy at the short waves is usually at a minimum, and vice versa. This 
suggests that energy exchanges take place continuously between the forced long waves 
and the transient disturbances. The disturbances remove energy from the planetary 
scale during the growing stage and feed energy back into the planetary scale during the 
decaying stage. As a consequence the amplitude of the time mean planetary-scale 
waves is smaller than that produced by a linear model. Since no dissipation is provided 
for the transient waves all energy extracted from the forced motion must either 
remain in the transient waves or be transferred back to the forced motion. Colton 
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(1973) speculated that the actual damping of the planetary scale by the local baro- 
tropic instability may be more effective if some dissipation mechanism for the 
transient synoptic waves is included. 

Although Colton's (1973) model lacks a sufficient damping mechanism and 
therefore produces a basic flow that is somewhat too intense and has a phase shift 
problem, it nevertheless illustrates the interesting phenomenon of waves developing 
and decaying in a spatially varying basic flow which contains locally unstable and 
stable regimes. To gain further physical insight into this phenomenon, especially with 
respect to the easterly jet of the summer monsoon, Tupaz, Williams and Chang 
(1977) formulated a 1 -level numerical model using the barotropic vorticity equation 
to study this problem from a linear instability point of view. In their model the basic 
zonal wind is approximated by a steady easterly Bickley (hyperbolic-secant-squared) 
jet which increases slowly from the inflow (eastern) boundary to a longitude in the 
middle of a rectangular domain, then decreases slowly downstream tov/ards the 
outflow (western) boundary. The basic meridional wind is derived in such a way that 
the total basic flow is non-divergent. Periodic forcing is introduced on the inflow 
boundary and a radiation condition, whereby a wave is allowed to propagate its 
energy out of the domain, is applied at the outflow boundary. The model is integrated 
from zero interior perturbation streamfunction until the perturbation fields vary 
periodically in time at each grid point. At this quasi-steady state the waves move 
through the region and grow’ or decay in space only, in reaction to the local stability 
properties of the basic flow. 

To compare the numerical results with the 'local grow’th rate' concept of the 
parallel flow theory, Tupaz et ciL (1977) developed a simple analytical model to 
calculate the effects of the local instabilities. This is done by considering an equation 
that allows for propagation and growth or decay, ij/[ -f c{x)\l/'^ = where / is 

time, X the zonal coordinate, and \p' the perturbation streamfunction. The local 
phase speed c and local growth rate // can be obtained from the parallel flow theory. 
When a periodic solution with a frequency (d and an amplitude F{x) is introduced, the 
downstream variation of F with respect to a reference point Xq may be written as 

Fix) = F{Xo) exp [-M dx/c(x)'] exp [ ^i(.y) ^/.v/V(a*)]. 

J xo J xn 

In this equation the local wavelength is o)/c and the local spatial growth rate is n/c. 
Figure 1 1 compares the steady state envelope for i// ' of the numerical solution and the 
corresponding F(x) for a case in which the central speed of the easterly jet is increased 
from 15 m sec' ^ at X = — 10 000 km to 30 m see' ^ at x = 0 and then reduced to 
15 m sec" ^ at x = 10 000 km. The curve for F(x) can be shifted up or down in the 
diagram by changing x^ , but its shape does not change. Two F{x) curves are included 
in Fig. 1 1, one for Xq = 280 km and the other for = 4200 km, both based on the 
local wavelength distribution measured from the numerical solution. It is evident 
from this figure that both the numerical solution and the local growth rate of the 
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X (1000km) 

Figure 1 1 

C omparison of wave packet envelopes in a locally unstable mean flow; (li/ computed from a numerical 
model ;F^(xo = 2S0 kmiandF^tv,^ = 4200 km), computed from an analytical model based on the parallel 

Mow theory. (From Tltaz et al. 19^7.) 



parallel flow theory put the maximum amplitude of the perturbations at .v = — 9800 
km. well past the jet maximum. This is in good agreement with the results obtained by 
Colton (1973) and Chang and Pentimonti (1977), and is qualitatively consistent 
with Krishnamurti and Roc'.ers’ ( 1970) observation that the most active disturbance 
area extends downstream from the jet maximum near the southern tip of India to 
central Africa (see Fig. 10.2, top panel, of Krishnamurti. 1971b). This longitude of 
maximum amplitude corresponds to the neutral point of the basic flow, and is 
immediately understandable from the local growth rate concept: a wave coming in 
from the inflow boundary will grow as it moves downstream in the locally unstable 
region, until it reaches the neutral point where it begins to decay. However, the 
numerical solution has a larger maximum amplitude than that predicted by the local 
growth rate of the parallel flow theory. The waves actually grow faster in the unstable 
region and decay faster in the stable region. Tupaz et al. (1977) evaluated the energy 
equation for the numerical model and found that the most important term for 
perturbation energy production is the Reynold stress term proportional to the 
meridional shear of the basic zonal wind. This term is also the only source term for a 
parallel flow model, and it depends on the phase structure of the perturbation field. 
Figure 12 compares the local phase tilt of the wave from the numerical model with 
that from the parallel flow theory at two longitudes: v = 0 km (Jet maximum) and 
A = —2800 km (slightly downstream from Jet maximum). At both places the phase 
tilt from the numerical model is significantly larger than that of the parallel flow^ theory. 
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Thus apparently the downstream variation of the basic flow augments the phase tilt 
of the waves which gives a larger growth rate in the unstable region. A similar situation 
is found in the stable region where the dynamic damping due to the continuous 
spectrum effect is enhanced by a stronger tilt in the same direction of the basic shear of 
the easterly jet. These results indicate that the barotropic energy exchanges associated 
with a zonally varying easterly jet which is locally unstable can actually be more 
prominent than those predicted by the local growth rate from the parallel flow theory. 
Hence the presence of the planetary-scale asymmetries during the summer monsoon 
strengthens the nonlinear energy conversions resulting in significant activity of 
synoptic-scale disturbances. 



Phase Tilt Angle 0 
-80 -60 -40 -20 0 20 



Phase Tilt Angle 0 ^ 





Figure 12 

Comparison of the phase tilts of waves in a locally unstable mean flow computed by a numerical model 
(^*) and the parallel flow' theory at two longitudes: (a) .x = 0 where r^lx, 0) = -30 m sec" ‘ and 
(b) .X = 2800 km where u(.x, 0) = -27.3 m sec" t From Tupaz t'/ i//. 1977.) 



Although Tupaz et al. (1977) employed linear equations only, their results may 
also shed some light on the effects of the transient waves on the basic flow. The waves 
remove kinetic energy from the basic flow (both zonal and planetary scale) and most 
of this energy is removed on the downwind side of the jet maximum. Therefore it may 
be expected that in a nonlinear extension of their model the planetary scale will 
experience maximum damping west of the easterly jet maximum, which undoubtedly 
will affect the amplitude and phase of the forced planetary-scale waves. 
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4. Oscillations of the planetary-scale monsoons 

The summer monsoon circulation is known to exhibit many fluctuations in a 
variety of time scales within each season. The most notable fluctuation is the alter- 
nating active and break periods of the monsoon. During the break periods most of 
the circulation characteristics may diminish or even reverse over the broad monsoon 
region. Less dramatic variations occur more frequently than the active-break cycles, 
although sometimes they are limited to local regions. Figure 13 shows the power 
spectrum of the 1967 summer 200 mb streamfunction over a region of the Tibetan 
high, analyzed by Krishnamurti et ai (1973a). A pronounced oscillation periodicity 
is noted around 10-13 days and a secondary spectral peak occurs near 3 days. A 
similar spectral distribution for the tropospheric wind field of the 1962 summer was 
found by M. Murakami (1976), where the main peaks occur at 10-15 days and 4-5 
days. Both of these investigators felt that the shorter period oscillation may be a 
manifestation of the synoptic-scale w'estward propagating monsoon depressions. The 
longer period oscillation, on the other hand, appears to be related to the break monsoon 
and other broad-scale variations of the monsoon circulation. The 10-15 day periodi- 
city of this oscillation has also been noted during many other years. For example, it is 
reflected in the daily rainfall data of 1963. 1971 (Murakami, 1972), and 1972, 1973 
(Krishnamurti and Bhalme, 1976). 

After detecting a spectral peak near the 13-day periodicity in many parameters of 
the planetary-scale monsoon system of the 1967 summer. Krishnamurti and Bhalme 
(1976) calculated the phase relationship between the oscillation of these parameters. 
Their result suggests that in a 13-day cycle the sequential strengthening of the mon- 




Figure 13 

Percent variance spectrum of the intensity of the Tibetan High during June August 1971. The dashed line 
IS a smoothed envelope. (From Krishnamurti ct ai.. 1973a.) 
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soon system including the monsoon trough, the cloud brightness, the surface Mas- 
carene high, the upper Tibetan high and easterly jet, and the central Indian rainfall 
occur after the dry and moist destabilization of the lower troposphere. This will then 
be followed by a weakening of the system in the same sequence after dry and moist 
stabilization. On the other hand, M. Murakami (1976) used a composite technique 
to study the 1962 summer monsoon, noting that it is difficult to verify the statistical 
significance of the cross-spectral quantities of this period range derived from a 
3-month time series of daily data. He found that the intensification and weakening of 
convection and other low- and upper-level winds are closely in phase for the 15-day 
cycle during this season. In any case, there is strong evidence that many components of 
the planetary-scale monsoon are involved in this oscillation. 

The observed 10-15 day oscillation raises an interesting question concerning its 
physical mechanism. Krishnamurti and Bhalme (1976) proposed that it may be a 
result of a cloud-radiative feedback mechanism in which the enhanced convection 
causes increased cloud cover which in turn reduces the incoming short-wave radiation 
and stabilizes the lower troposphere. This is then followed by decreased convection 
and cloudiness, and therefore increased short-wave radiative heating at the surface 
which renews the active convection. Another possible oscillation mechanism is 
revealed by a simple, domain-averaged, atmosphere-ocean interacting model by 
Webster and Lau ( 1977). In their model the large heat capacity of the ocean causes 
it to respond slowly to the solar heating while the atmosphere temperature reaches a 
maximum due to heating at the land surface. After a period of time the ocean tem- 
perature peaks and this feeds back to the atmosphere through convection resulting in 
another temperature maximum. Although this process alone is unlikely to produce 
subsequent maxima without the participation of other mechanisms, it may still be 
relevant to some of the active-break cycles because Rama.murthv's (1969) survey 
(quoted in Krishnamurti and Bhalme, 1976) indicates that there are only one or two 
major breaks during most of the summer monsoon seasons. Other possible mechan- 
isms for the planetary-scale oscillation may include interactions between hemispheres 
and between stationary and propagating long waves as mentioned by Krishnamurti 
and Bhalme (1976), and nonlinear vacillations between different scales. The latter 
mechanism produces a 5-day oscillation of the planetary scale in Colton's (1973) 
barotropic model, which is purely dynamical and does not provide dissipation 
processes for the synoptic scale. If thermodynamics and synoptic damping are 
included, the interaction would become more complex and the oscillation period 
could be modified. 

The planetary-scale circulation during the winter monsoon also exhibits significant 
variations, largely due to the influence of the cold air outbreaks from the Asian 
continent. These outbreaks often interfere with the main winter convective region in 
the maritime continent and cause widespread weather changes in the South China 
Sea and Malaysia-lndonesia region, including the development of strong synoptic- 
scale disturbance activity. Since this convective region provides the major energy 
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source for the thermally direct, planetary-scale east-west overturning and Hadley 
cells (Krishnamurti ei al. 1973b), its variations directly influence the circulations in 
the entire tropics and southern mid-latitudes during winter. These planetary-scale 
zonal and lateral teleconnections are reflected in the analyses by Murakami and 
Unninayar (1977) and Krishnamurti et al. (1973b) despite their limited data base. 
The initiation and progress of theoretical studies in this area would depend heavily on 
the expansion of our present observational knowledge of this problem. 



5. Concluding remarks 

The three problems reviewed here are inter-related to each other. For example, 
the localized barotropic instability is a reflection of the strong nonlinear energy 
transfer which appears to be quite relevant to the damping of the planetary-scale 
flow during the summer monsoon. Also Chang and Pentimonti's ( 1977) experiment 
indicates that these energy conversions are more pronounced when the planetary 
scale contains an appreciable oscillation. The enhanced upper tropospheric in- 
stability in summer due to both time and spatial variations of the planetary scale also 
has important ramifications for the synoptic-scale disturbances. At least in the mid- 
occanic trough regions, the upper-level disturbances may develop in the vertical and 
could be related to some of the weather-producing easterly waves and tropical 
cyclones. A different type of problem arises for the winter monsoon during which the 
upper tropospheric synoptic-scale waves and zonal flow feed energy to the planetary 
scale (Krishnamurti ct al. 1973b). This on the one hand removes the nonlinear 
processes as a possible damping mechanism, and on the other hand requires an 
explanation for the energy source of the synoptic scale. Clearly much more effort is 
needed to answer these questions, whose importance lies in the fact that a thorough 
understanding of the planetary scale is crucial for the proper simulation of the entire 
multiple-scale monsoon circulation. 
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